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Abstract 

We reexamine the issue of the soliton mass in two-dimensional models with M = 1 
supersymmetry. The superalgebra has a central extension, and at the classical level the 
soliton solution preserves 1/2 of supersymmetry which is equivalent to BPS saturation. 
We prove that the property of BPS saturation, i.e. the equality of the soliton mass 
to the central charge, remains intact at the quantum level in all orders of the weak 
coupling expansion. Our key finding is an anomaly in the expression for the central 
charge. The classical central charge, equal to the jump of the superpotential, is amended 
by an anomalous term proportional to the second derivative of the superpotential. The 
anomaly is established by various methods in explicit one-loop calculations. We argue 
that this one-loop result is not affected by higher orders. We discuss in detail how the 
impact of the boundary conditions can be untangled from the soliton mass calculation. 
In particular, the soliton profile and the energy distribution are found at one loop. A 
"supersymmetry" in the soliton mass calculations in the non-supersymmetric models is 
observed. 
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1 Introduction 



The issue of quantum corrections to the soliton mass in two-dimensional models with 
M = 1 supersymmetry has a long and dramatic history. As was noted in Ref . |lj , in the 
models with central extensions topologically stable solitons can be Bogomol'nyi-Prasd- 
Sommerfeld (BPS) saturated. Such solitons are annihilated by the action of 1/2 of the 
supergenerators, so that 1/2 of the supersymmetry (SUSY) is preserved. The remaining 
1/2 of the supergenerators act non-trivially. The mass of the BPS saturated solitons 
must be equal to the central charge. The simplest two-dimensional model with the min- 
imal superalgebra, admitting solitons, was considered in 0. This is the Wess-Zumino 
model with one real scalar field and one two-component real spinor (to be referred to as 
the minimal model hereafter) . It was argued that, due to a residual supersymmetry, 
the mass of the soliton calculated at the classical level remains intact at the one-loop 
level. A few years later it was noted || that the non-renormalization theorem 0] can- 
not possibly be correct, since the classical soliton mass is proportional to m 3 /A 2 (where 
m and A are the bare mass parameter and coupling constant, respectively), and the 
physical mass of the scalar field gets a logarithmically infinite renormalization. Since 
the soliton mass is an observable physical parameter, it must stay finite in the limit 
M uv — > oo, where M uv is the ultraviolet cut off. This implies, in turn, that the quantum 
corrections cannot vanish - they "dress" m in the classical expression, converting the 
bare mass parameter into the renormalized one. The one-loop renormalization of the 
soliton mass was first calculated in Technically the emergence of the one- loop cor- 
rection was attributed to a "difference in the density of states in continuum in the boson 
and fermion operators in the soliton background field" . The subsequent work Q dealt 
with the renormalization of the central charge, with the conclusion that the central 
charge is renormalized in just the same way as the soliton mass, so that the saturation 
condition is not violated. 

Since then a number of one-loop calculations have been carried out |5|-|T5|. The 



results reported and the conclusion of saturation/non-saturation oscillate with time, 
with no sign of convergence. Although all authors clearly agree that the logarithmically 
divergent term corresponds to the renormalization of m, the finite term comes out 
differently, varying from work to work, sometimes even in the successive works of the 
same authors, e.g. |13|,|l4l or |||9|]. According to a recent publication |fTI| , the BPS 



saturation is violated at one loop. This assertion reverses the earlier trend P, |T0|.|TT 



according to which the soliton in the minimal model is BPS saturated: the central 
charge and the soliton mass are renormalized in a concerted way. 

The problem becomes obscure partly due to the fact that some authors (e.g. ||) 
use non-supersymmetric boundary conditions. When non-supersymmetric boundary 
conditions are imposed, explicit supersymmetry is lost. Gone with it is the cancellation 
of the boson and fermion loops inherent to supersymmetry. Then one has to consider 
separately the boson and fermion contributions - each of them individually diverges. 
Extracting finite terms in the difference of two diverging contributions complicates the 
analysis. This approach - the separate treatment of the boson and fermion corrections 
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in a SUSY-insensitive manner - prevails in the recent works. Although such approach 
is fully legitimate, it does not take advantage of SUSY. 

We reexamine the issue of the quantum corrections to solitons in two-dimensional 
models with minimal supersymmetry. Our strategy is to take maximal advantage of 
super symmetry, at every step. One of our main goals is establishing a relation with other 
similar problems where 1/2 of SUSY is preserved in the background field: the instanton 
calculations in SUSY gluodynamics |l6j and the BPS-saturated domain walls |17|,^8|. 
In both cases the quantum corrections totally cancel, bosons against fermions. Up to a 
point, a similar cancellation takes place in the two-dimensional models too. If it were not 
for boundary effects, the classical expression for the soliton mass would not be corrected 
at the quantum level. The boundary effects, which play no role in SUSY instantons or 
the BPS-saturated domain walls, are responsible for the soliton mass renormalization in 
the two-dimensional models with minimal supersymmetry. The divergence between the 
two- and four- dimensional problems lies here. A non-renormalization theorem survives 
in the form of the statement that the soliton mass is equal to the central charge, to 
all orders. We present a general argument for why the BPS saturation cannot be 
violated by quantum corrections in the weak coupling regime. In four dimensions the 
phenomenon known as multiplet shortening protects the BPS states from becoming 
nonsaturated at higher orders. In two dimensions the multiplet shortening does not 
occur, literally speaking. However, the multiplet shortening per se is not the reason 
lying behind the protected nature of the BPS states, but, rather, a signature. The 
genuine reason is the existence of an invariant subalgebra of the superalgebra, i.e. the 
existence of a supercharge Q possessing the following properties: (i) Q 2 is nontrivial; 
(ii) Q annihilates all states from the given multiplet. Although the multiplet shortening 
does entail the existence of such a subalgebra, the inverse is not necessarily true. In 
fact, the two-dimensional models with minimal supersymmetry present an example of 
the problem where the invariant subalgebra of the required type exists without the 
multiplet shortening. This is sufficient for establishing the theorem "soliton mass = 
central charge" , order by order. 

The next question to address is: "what is the central charge in the theories with 
minimal supersymmetry?" This is a key element of our analysis. As we will show below, 
much of the confusion in the current literature is due to a tacit assumption that the 
central charge Z is equal to the difference between the values of the superpotential W 
at spatial infinities, 

Z = W[<j>(z = oo)] - W[(f)(z = -oo)] . (1.1) 

This expression was obtained long ago M . At the classical level it is a valid expression. 
It was implied in the previous studies that the only change at the quantum level is the 
substitution of the classical superpotential W in Eq. (|1.1| ) by its operator form. We 
show that the central charge is additionally modified by a quantum anomaly, 

w 

W^W + — , (1.2) 
i.e. the anomalous term W" should be added in the expression for the central charge. 
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We have obtained this anomaly by four different methods: (i) using dimensional 
regularization; (ii) regularization by higher derivatives; (iii) from an effective super- 
potential; (iv) from the models with the extended (Af = 2) supersymmetry deformed 
by soft terms breaking of M = 2 down to M = 1. If the first two methods are quite 
straightforward, the third and the fourth require comment. As in many other similar 
problems, it is convenient to calculate the soliton mass starting from an effective La- 
grangian (superpotential), including the loop effects by construction, rather than from 
the bare Lagrangian. It is not difficult to obtain the effective superpotential at one-loop 
level, which results, in turn, in a certain prediction for the soliton mass. This predic- 
tion can be compared with the result obtained through the central charge. The two 
calculations agree with each other provided one accounts for the anomaly in Eq. (|1.2j ). 

An alternative approach is based on the models with extended supersymmetry, 
M = 2. The M = 2 theories in two dimensions can be obtained by dimensional reduction 
of the Wess-Zumino model , from four to two dimensions. The superpotential is not 
renormalized [2(J, and so neither is the soliton mass. The boundary effects - the prime 
culprits for the soliton mass renormalization in the minimal model - vanish in the 
M = 2 case. We then deform the M = 2 model by introducing terms which softly break 
M = 2 SUSY down to M = 1. The emerging dynamics is quite rich and interesting. On 
the technical side, the softly broken M = 2 theories provide a possible regularization 
procedure for M = 1 theories. Applying this procedure we recover Eq. (|1.2p again, for 
the fourth time. 

Although the fact that the central charge does have anomalies is very well known in 
the four-dimensional theories (e.g. the existence of the anomaly in two-dimensional 
theories with minimal supersymmetry was overlooked in all previous publications. 
Needless to say that the additional term proportional to W" in Eq. ( |1 . 2|) gives a 
non-vanishing contribution to the central charge. In conjunction with the statement 
of BPS saturation, this contribution translates into a (finite) extra term in the soliton 
mass. If the BPS saturation per se is established on the basis of comparison between 
two distinct calculations - of the soliton mass, on the one hand, and of the central 
charge on the other - the omission of the anomalous term leads to an apparent para- 
dox. The general theorem forbidding the loss of BPS saturation in quantum loops in 
the weak coupling expansion seems to be violated. 

Since the quantum corrections are associated with boundary effects, their determina- 
tion reduces to a calculation in the "empty" vacuum, which can be trivially performed. 
At the same time, the standard calculation of the soliton mass becomes sensitive to 
the choice of boundary conditions. This is another subtle point which was not properly 
treated in the literature. It was believed that the soliton mass depends on the boundary 
conditions, and that particular choices of the boundary conditions are "more correct" 
than others. It is clear that physically the soliton mass is saturated in the domain near 
its center, and must not depend on the boundary conditions set at a large distance from 
the soliton center. 

Although the boundary conditions do not affect the physical soliton mass, some 
choices are practically more convenient. In fact, inappropriate choices contaminate the 
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field near the boundary and may shift the total energy in the box by a quantity which is 
of the same order as the "genuine" one-loop quantum correction to the kink mass. If the 
boundary domain is not treated properly, the resulting calculation yields the mass of a 
contaminated soliton, rather than the physical mass of the ground state in the sector 
with unit topological charge (i.e. the true soliton). Much of the controversy existing 
in the literature is due to this contamination. We explicitly demonstrate how these 
distortions appear in the boundary domain, and show how the boundary conditions 
can be untangled from the soliton mass calculations. One of the possibilities to control 
the contamination is to consider the energy- density distribution in the soliton at one 
loop. This is done in Sec. [5| In addition, we perform a calculation of the one-loop 
correction to the soliton profile. 

The last issue addressed in this work is the soliton mass calculation in the non- 
super symmetric sine-Gordon model. Surprisingly, this problem is naturally formulated 
in super symmetric terms too, although the original Lagrangian has no explicit super- 
symmetry. The role of the SUSY partners is played in the case at hand by the sectors 
with the zero and unit topological charges, respectively. This approach yields a num- 
ber of useful results. In particular, closed analytical expressions for the boson and 
fermion Green functions in the soliton background, known previously |21,Effi|, are given 
a transparent interpretation. 

The outline of the paper is as follows. In Sec. [| we briefly review the formulation 
of the problem, and general aspects of two-dimensional theories with minimal super- 
symmetry. Most of the results presented here are not new. Since they are scattered in 
the literature we find it convenient to collect them in one place, organizing a common 
perspective. We briefly review the superfield formalism in two dimensions and describe 
the two models with minimal supersymmetry we will deal with: the polynomial model 
and the sine-Gordon model. Then the kink solutions are discussed, and the notion of 
the kink superfield is introduced. A residual supersymmetry which annihilates the kink 
is identified. At the classical level, the corresponding SUSY generator acts on the kink 
trivially. It is shown that this property cannot be lost in perturbation theory. In Sec. 
[2.3| we discuss Witten's index in two-dimensional models with minimal supersymmetry. 

Section ^|is devoted to derivations of the central charge. The focus is on the anoma- 
lous part of the central charge, which is obtained in a variety of ways. We formulate a 
convenient (and universal) regularization procedure based on higher derivatives, which 
regularizes all ultraviolet divergences. We also explain how the same anomaly emerges 
as an infrared effect, through the level flow. This is a standard situation - the anomaly 
has two faces: ultraviolet and infrared (see e.g. [f23fl). In Sec. a local version of the 
BPS saturation condition is explored. In Sec. ^ the one-loop corrections to the kink 
mass are calculated, both in general form and in the two models under consideration, 
the polynomial model and sine-Gordon. We present two distinct calculations. Section 
[| deals with the soliton profile at one loop. The energy- density distribution in the kink 
is also presented here. The results are based on the explicit expressions for Green's 
function in the soliton background. In Sec. |6] the calculations of the soliton masses 
in non-supersymmetric polynomial and sine-Gordon models are reformulated in terms 
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of a supersymmetry which - surprisingly - can be introduced in these problems. Sec- 
tion presents the calculation of the soliton mass in the supersymmetric sine-Gordon 
model at two loops. Finally, Sec. |8| is devoted to extended supersymmetry and Af = 1 
deformations of the M = 2 model. Section [| summarizes our findings. 

2 The minimal SUSY in two dimensions 

2.1 Superspace for J\f=l SUSY models 

The two-dimensional space x^ = {t, z} can be promoted to superspace by adding a 
two-component real Grassmann variable 9 a = {9i, 82}- The coordinate transformations 

a ^9 a + e a , x» -> x^ -iO^e (2.1) 

add SUSY to the translational and Lorentz invariances. A convenient representation 
for the two-dimensional Majorana 7 matrices we use throughout this paper is 

7° = °2 , 7 1 = 1V3 • (2.2) 
The real superfield $ (r, 9) has the form 

$ (x, 9) = 4) (x) + 9^{x) + -99F(x) , (2.3) 

where 9 = #7°, and 9, ip are real two-component spinors. The superspace transforma- 
tions ( |2.1| ) generate the following SUSY transformations: 

5<j) = sip , 5%l) = -i <9 M Ye + Fe, 5F = -ie^d^ . (2.4) 

In 1+1 dimensions it is a minimal, M = 1, supersymmetry. 

The action of a model invariant under the transformations (|2.4|) is 

S = ijd 2 9 d 2 x i-D a <S>D a <& + W($)| (2.5) 

where W($) will be referred to as the superpotential, keeping in mind a parallel with the 
four-dimensional Wess-Zumino model, although in the case at hand the superpotential 
term is the integral over the full superspace, and is not chiral. Moreover, d 2 x = dt dz; 
the spinorial derivatives are defined as follows 

D a = w -i (Y9) a d, , D a = -— + % , (2.6) 

so that 

{D a D p } = 1% tf^apd^ . 
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In components the Lagrangian takes the form 



£ = o d ^ + + F 2 ) + W'(0)F - -W"((f))il>il> . (2.7) 

In two dimensions any superpotential function W leads to a renormalizable field 
theory and is thus allowed. Many of the results presented below are independent of 
the particular choice of W. For our purposes it is generally sufficient to limit ourselves 
to odd functions of i.e. W($) = — W(— $) (Sec. |] is exception). We require the 
superpotential W(0) to have more than one extremum, i.e. the equation dW/d0 = 
must have at least two solutions. Then one deals with several degenerate vacua. 
The field configurations interpolating between distinct vacua at the spatial infinities 
(z — > ±oo) are topologically stable solitons (kinks). 

In certain instances we find it illustrative to consider specific examples. Two choices 
are of most practical interest: the polynomial superpotential, and the minimal super- 
generalization of the sine-Gordon model. In the first case the superpotential function 
W is parameterized as follows: 

W($) = — $ - £ $ 3 , (2.8) 
where the parameters m and A are real numbers. In the second case 

W($) = mv 2 sin - . (2.9) 

v 

The first model will be referred to as SPM, the second as SSG. 



2.2 Multiple vacua and kinks 

The SPM model has two vacuum states corresponding to <ftf = ±ra/(2A). A kink 

solution interpolates between 4>~ = — m/2A at z — >• — oo and 0+ = m/2A at z — > oo, 
while an antikink between 0+ = m/2A and <p~ = —m/2X. The classical kink solution 
has the form 

O = — tanh— . (2.10) 
In the sine-Gordon model there are infinitely many vacua; they lie at 




(2.11) 



where k is an integer, either positive or negative. Correspondingly, there exist solitons 
connecting any pair of vacua. In this case we will limit ourselves to consideration of 
the "elementary" solitons connecting the adjacent vacua, e.g. = ±nv/2, 

0o = f arcsin [tanh(mz)] . (2-12) 
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Both solutions ( |2.10| ) and ( |2.12| ) are BPS saturated, i.e. they satisfy the first order 
differential equation, 

**M = -F{4> ) = W{<l>a). (2.13) 

In terms of superfields the BPS saturation means that a background field $o is of the 
form 
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(2.14) 



where / is some function. In other words, instead of the generic superfield <3>(x, 9) 
depending on two coordinates x M and two Grassmann coordinates 9 a we deal with 
a background (super)field depending only on one specific combination, z — 88/2. In 
components Eq. ( |2.14j ) implies that 



4>o = f(z), 



df(z) 
dz 



(2.15) 



The soliton solutions ( |2.10|) and (|2.12|) belong to this class, with a special choice of the 
function f(z). 

A special, topological nature of the background field ( 2.14 ) is seen from the following 
property: the total energy is defined by the asymptotics at z — > ±oo, 

- j £[%}dz = AW = W[<j)(z = oo)] - W[<p(z = -oo)] , (2.16) 

irrespective of the particular choice of the function /. Substituting the asymptotic 
values of <fi at ±oo in Eq. ( |2.16| ) we find the soliton mass 



m 
6 A 2 

2 mv 2 



(SPM) 
(SSG) 



(2.17) 



where the subscript emphasizes that this result is obtained at the classical level. 

For any background field of the type ( J2.14 ) half of the supersymmetry is preserved. 
Indeed, let us consider supertransformations (|2.4| ) with the parameter e<i = 0, 

(2.18) 



S8i = ei , 58 2 = , 5t = —i8\Ei , 5z = i8 2 ei . 
This transformation leaves the interval 

z inv = z — \-88 = z + idid 2 



(2.19) 



invariant. Consequently, backgrounds of the form ( |2.14|) are not changed. It means 
that quantization preserves the invariance under the following transformations of the 
fields, 
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"#2 



E X 



-icb 



d z d> + F 



E X 



-iipx 



(2.20) 
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This form shows the multiplet structure of the surviving supersymmetry. An important 
consequence of the invariance ( |2.20|) is that the fermion-boson degeneracy is not lost for 
quantum excitations although a part of SUSY is broken in the given fixed background. 

The part of supersymmetry which is spontaneously broken by the soliton configu- 
ration corresponds to the transformations produced by 62, 

50! = , 59 2 = e 2 , 5t = -i9 2 S2, 5z = i6 1 s 2 . (2.21) 



Applied to the soliton solution it generates the fermionic field ip2 (use Eq. (|2.4j ) with 

si = oy. 

6^2 = (-d z <t> + F ) e 2 = -2 e 2 • (2.22) 

dz 

In addition to the above fermion zero mode there is also the bosonic zero mode generated 
by the spatial translation. Both zero modes can be conveniently described by the 
introduction of two collective coordinates: the soliton center zq and the real Grassmann 
parameter 77 corresponding to the fermion zero mode, in much the same way this was 
done for instantons and domain walls [ I8| . To this end let us make the substitution 

z imi = z + 16162 — > z- mv = z - z Q + 161(62 - rj) . (2.23) 

The kink superfield, with all collective coordinates included, is 

$o = /(*w), (2.24) 

where the function / is given by Eq. (|2.10|) for the polynomial model and by Eq. Q2.12 ) 
for the sine-Gordon one. 

The deformation of the kink superfield $0 under the supertransformations ( |2.1|) is 
compensated by the following transformations of the collective coordinates: 

5z = ir] ex, 5r] = 26 2 - (2.25) 

The parameters z and 77 can be treated as time-dependent dynamical variables, 
describing the superkink. The quantum mechanics of these coordinates emerges upon 
integration over all non-zero modes in the field-theoretical model at hand, 



V[(/)(x)]V[tp(x)] exp{z / £(0, tp)d 2 x} 



V[z (t)]V[rj(t)] exp \i dt 



-M + -Mz\ + —7777 



(2.26) 



where the functional integral is performed in the sector with unit kink number, and M 
is the physical mass of the kink. 
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2.3 Witten's index in the models with minimal supersymme- 
try 

Here we comment on a general aspect of the two-dimensional models with minimal 
supersymmetry which is only indirectly related to the solitons. In the model (|2.5| ), 
( f2.7| ) with the polynomial potential the Witten index is zero. This can be seen in many 
ways. Say, if the mass parameter m 2 in the superpotential (|2.8| ) is made negative, the 
equation dW/d$ = has no real solutions. At negative m 2 the Z2 symmetry 

$->-<!>, i> -» 75^ 

stays unbroken. This protects <j) from developing a vacuum expectation value. There- 
fore, the ip held is massless, and is the Goldstino of the spontaneously broken SUSY. A 
convenient order parameter is (W), where 

W'^ = ^-A$ 2 . (2.27) 
d$ 4A 

Since W = (1/2)DD<&, the expectation value (W) = in the supersymmetric vacuum. 

Alternatively, one can consider the theory with positive m 2 in a finite box. If, 
following Witten ||25|| , we make the box width L small, XL <C 1, and retain only the 



zero momentum modes, discarding all others, we get a quantum-mechanical system 
known to have a vanishing Witten index [EH . 



At positive (and large) values of m 2 the Z2 symmetry is spontaneously broken, so 
that the theory has two physically equivalent vacua. The vanishing of the Witten index 
implies that if one of the vacua is bosonic, the other is fermionic. This is a unique feature 
of the minimal model: in four-dimensional field theories and in two-dimensional field 
theories with extended supersymmetry (and no massless particles) all vacuum states 
are clearly of the bosonic type. 

All theories with vanishing Witten index are potentially unstable with respect to 
spontaneous SUSY breaking. Moreover, in a finite volume this spontaneous break- 
ing does occur. However, in the infinite volume limit spontaneous SUSY breaking is 
impossible provided X 2 /m 2 <C 1. Indeed, one can show that X 2 /m 2 is the genuine di- 
mensionless expansion parameter in the theory at hand. At small X 2 /m 2 the theory is 
weakly coupled. Since there are no massless fields in the Lagrangian, (and they cannot 
appear as bound states in the weak coupling regime), there is no appropriate candidate 
to play the role of Goldstino, and, hence, supersymmetry must be realized linearly. 

The Witten index for the sine-Gordon model also vanishes. This can be checked 
by putting the theory in a finite box. The resulting quantum mechanics has two su- 
persymmetric vacua - one bosonic and one fermionic - provided the wave functions 
considered are periodic. One can lift these two vacua from zero by imposing Bloch 
boundary conditions on the wave functions, with the angle 9 7^ 0. 
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3 Ultraviolet and infrared regularization 

3.1 Superalgebra (unregularized) 

We start with the operator construction which corresponds to the classical level of 
consideration. The SUSY algebra is built on the supercharges Q a , defined as 

Q a = JdzJ°, (3.1) 

where the conserved supercurrent is 

= (^0)7 V0 - i-FW . (3.2) 
Using the canonical commutation relations [] 

}(t, z),<j>(t, z')]=i6(z-z'), {Mt, z),Mt, z')} = ^)^5{z-z') (3.3) 

one finds the anticommutator, 

{ J£, Qp) = 2 { lv ) aP + 2i ( 7 5 ) a/3 ^ , (3.4) 
where is the energy-momentum tensor 

^ = d u 4> + -fyf idty - - [<9 7 <9> - F 2 ] , (3.5) 

and C is the conserved topological current, 

^ = e ^d u W . (3.6) 

Symmetrization (antisymmetrization) over the bosonic (fermionic) operators in the 
products is implied in the above expressions. Also, the notation, 7 5 = 7°7 1 = — <7i 
and F = — W, is used. 

Integrating the \i = component of Eq. ( |3.4| ) over space gives the SUSY algebra: 

{Q a , Qfs} = 2 ( 7 M ) a/3 P M + 2t ( 7 5 ) a/3 Z . (3.7) 

Here P M = /dz^o^ are operators of the total energy and momentum, and Z is the 
central charge, 

Z = Jdz(° = Jdz d z W((f)) = W[<f>(z = oo)] - W[(j)(z = -oo)] , (3.8) 

which coincides with the topological one. In the theories ( |2.7| ) with several degenerate 
vacua the central charge Z is nontrivial for the soliton configurations. Thus, the M = 1 
superalgebra gets centrally extended [flj. 

Let us emphasize that Eq. ( |3.4|) establishes supersymmetry relations between the su- 
percurrent J M , the energy-momentum tensor fi^, and the topological current and for 
this reason cannot be changed by quantum corrections. The explicit expressions ( |3.2| ), 
( ^5|) and (|3.6j ) for i?^ and are classical ones, quantum corrections will change them. 



^^Here and below <f> should be understood as the canonical momentum ir^ conjugate to 
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3.2 Anomaly in the central charge 

The structure of the centrally extended superalgebra ( |3.7| ) is quite general. At the same 
time, the particular expression for the central charge depends on the dynamics of the 
model under consideration. Equation ( |3.8| ) is obtained at the classical level. It is easy 
to see that quantum corrections do modify it. Here we will present a simple argument 
demonstrating the emergence of an anomalous term in the central charge. We also 
discuss its physical meaning. A detailed direct calculation will be carried out after we 
introduce the ultraviolet regularization by higher derivatives. 

To begin with, let us consider 7 /i J M where J M is the supercurrent defined in Eq. ( |3.2j ). 
This quantity is related to the superconformal properties of the model under consider- 
ation. At the classical level 

(7%) class = -2ii^. (3.9) 

Note that the first term in the supercurrent (|3.2j) gives no contribution in Eq. (|3.9| ) due 
to the fact that in two dimensions 7 M 7 l/ 7 /i = 0. 

Multiplying Eq. (|3^ ) by 7 M from the left we get the supertransformation of 7 M J M , 

\ {7%,<5}=^ + W 5 C M . (3-10) 

This equation establishes the supersymmetric relation between 7 M J M , ^ and and, 
as was mentioned above, remains valid with quantum corrections included. But ex- 
pressions for these operators could be changed. Classically the trace of the energy- 
momentum tensor is 

as follows from Eq. (|3.5| ). The zero component of in the second term classically 
coincides with the density of the central charge, d z W, see Eq. (|3~^). It is seen that the 
trace of the energy-momentum tensor and the density of the central charge appear in 
this relation together. 

It is well-known that in renormalizable theories with ultraviolet logarithmic diver- 
gences, both the trace of the energy- momentum tensor and 7 M J M have anomalies. We 
will use this fact, in conjunction with Eq. (|3.10|) , to establish the general form of the 



anomaly in the density of the central charge. 

To get some idea of the anomaly, it is convenient to use the dimensional regulariza- 
tion. If we assume that the number of dimensions D = 2 — e rather than D = 2 the 
first term in Eq. ( |3.2| ) does generate a nonvanishing contribution to 7 M J M , proportional 
to (D — 2)(d u (p) r f^). At the quantum level this operator gets an ultraviolet logarithm 
(i.e. {D — 2) _1 in dimensional regularization), so that D — 2 cancels, and we are left 
with an anomalous term in 7 M J M . 

To do the one-loop calculation we apply here (and in many instances below) the 
background field technique: we substitute the field by its background and quantum 
parts, (j) and x, respectively, 

— <P + X- (3.12) 



12 



Specifically, for the anomalous term in 7 M J^, 



(7 M ^)anom= ( D ~ 2 ) W>) 7*V = -{D-2)jcfd^ = i (D - 2) X W (0 + x) * , (3.13) 

where an integration by parts has been carried out, and a total derivative term is 
omitted (on dimensional grounds it vanishes in the limit D = 2). We also used the 
equation of motion for the ip field. The quantum field \ then forms a loop and we get 
for the anomaly, 

(7%) anom = i(D-2) (0| X 2 |0) W'"(0) * = -{D - 2)/^j W"'(0) V> 

= -^W'"(0)^ (3.14) 

Z7T 

The supertransformation of the anomalous term in is 

^ {(7%) anom >Q} = (^ - i^"' F ) + *rrf<™» (i^") • (3-15) 

The first term on the right-hand side is the anomaly in the trace of the energy- 
momentum tensor, the second term represents the anomaly in the topological current, 
the corrected current has the form 

^ = e ^d v (W + ^- W") . (3.16) 



An 

Consequently, the operator of the central charge 

Z=(w + ^-W") -(W + —W") , (3.17) 

is the difference of W + (W"/47r) at z = ±00. Other anomalies in the supermultiplet 
together with their classical counterparts have the following form, 

7% = 2i[W + —)il>, 

( W"\ 1 / W"\ - 

k = - F + 17J + 2 ( w " + ^rj ^ ■ < 3 - 18 ) 

The physical meaning of the anomalous term in the central charge becomes quite 
apparent from the consideration above. The supersymmetry relates the anomaly in the 
topological current with a response of the system to the scale transformations, 

x^(l-p)x, 9 -»• (1- {p/2))6, (3.19) 



The kinetic term of the action (|2.5|) is invariant. The non-invariance of the classical 
superpotential term is associated with W. This is not the end of the story, however. 
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To reveal the anomaly it is convenient, simultaneously with ( |3.19| ), to rescale the pa- 
rameters in the superpotential in accordance with their canonical dimensions. For 
instance, in the SPM model m — > (1 + p)m , A — (1 + p)A, and in the SSG model 
m — > (1 + p)m , v — > v. Then, 

W->(l + p)W, (3.20) 

and the classical action is invariant (the canonical dimension of W in two dimensions is 
unity). The quantum action is not invariant, due to the occurrence of a new dimensional 
parameter, the ultraviolet regulator mass, see e.g. Eq. ( |4.15|) . As it follows from Eq. 
( |4. 15| ) , the non- invariant contribution is W"/47r. This term gives rise to the anomaly 
in the central charge density. 

It is instructive to compare the above pattern with the anomaly in the central charge 
in the four-dimensional Wess-Zumino models, discussed previously in the literature [18 . 
In four dimensions, the central charge is proportional to the jump in 

- E + I E , (3-21) 

where the summation runs over all chiral superfields $j on which the given Wess-Zumino 
model is built. Here 7, is the anomalous dimension of the field $j. The first two terms 
are purely classical. The coefficient 3 in front of W is due to the fact that the canonical 
dimension of the superpotential in the four- dimensional Wess-Zumino model is three. 
The second term reflects the canonical dimension of the superfield $ which is equal to 
unity in four dimensions. (This is in contrast to two dimensions, where the canonical 
dimension of $ vanishes.) The first two terms in Eq. ( |3.21| ), combined together, vanish 
for purely cubic superpotentials. The reason is quite clear - if the superpotential in the 
four- dimensional Wess-Zumino model is a cubic polynomial, such a model is classically 
scale invariant, so the classical part in the central charge must also vanish. 

The last term in Eq. (|3.21|) is the quantum anomaly. It reflects the fact that the field 
$ in loops acquires an anomalous dimension, through a logarithmic renormalization of 
the kinetic term. In fact, due to the non-renormalization theorem |HJ, the Z factor 
for the kinetic term is the only renormalization occurring in the theory. Hence, the 
term proportional to the anomalous dimensions 7^ is the only anomalous term. In two 
dimensions, the kinetic term receives no logarithmically divergent renormalizations, 
and the anomalous dimension of $ vanishes. However, unlike four dimensions, the 
superpotential itself gets renormalized. This is the origin of the W'/Att anomalous 
term in Eq. ( |3.16|) in two dimensions. In four dimensions the anomalous term is a 
full superderivative, its expectation value over the supersymmetric vacuum vanishes 
automatically. Therefore, it plays no role in the computation of the central charge for 



the domain walls [fLq| . For the very same reason, there is no anomalous correction in 
the central charge in the two-dimensional soliton problems with M = 2. At the same 
time, the expectation value of the W" anomalous term in two dimensions by no means 
vanishes. The contribution of the anomaly to the central charge (and, hence, to the 
mass of the BPS saturated solitons) is absolutely essential for the overall self-consistency 
of the analysis. 
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3.3 Representations of the centrally extended algebra 



Let us consider the representation of the algebra (|3.7|) in the soliton sector. In the rest 
frame, where P z = 0, we have 

Ql = H + Z, Q\ = H-Z. (3.22) 

The supermultiplet consists of two degenerate states: the bosonic | sol;,) and fermionic 
| sol/) kinks.0 The supercharge Qi connects them, 

Qi|sol 6 ) = V2M|sol/) , Qi| sol/) = V2M| soU). (3.23) 

The second supercharge Q2 annihilates both soliton states, 

Q 2 I sol) = . (3.24) 

This equation is the operator form of the BPS saturation condition. Previously we 
discussed both features, Eq. ( |3.23| ) and Eq. ( |3.24| ), at the level of the classical solution. 

Note, that an analogue of the condition ( |3.24j ) in four dimensions is related to 
multiplet shortening. This multiplet shortening is possible only if there is a supercharge, 
which similarly to Q2, annihilates some states, as in Eq. (|3.24|) , and its square is related 
to H — Z, as in Eq. ( |3.22| ). In fact, the multiplet shortening is nothing but a signature 
(sufficient condition) of the existence of an invariant subalgebra. Although we do not 
have multiplet shortening in two dimensions, the invariant subalgebra formed by Q2 
and Q\ still exists. The relation ( |3.24j ) is enough to guarantee the absence of quantum 



corrections to (sol|i7 — Z \ sol). 

Can the soliton become non-saturated after the quantum corrections are switched 
on? The answer is no. In any supersymmetric theory, if Q2 annihilates a state at 
the tree level, it will continue to do so to any finite order if the coupling constant is 
small (for a review see e.g. p5]). The proof is based on the fact that the number of 
degenerate states in the multiplet cannot be increased from two to four in the weak 
coupling regime. Indeed, if both supercharges, Q\ and Q2, acted nontrivially one could 
construct two fermion states, Qi|solfe) and Q 2 |soU). Then the second bosonic state 
would be obtained as sol^). Since at the classical level there is no extra pair of 

degenerate states it cannot appear in perturbation theory. Thus, the soliton is protected 
from losing its saturated nature. 

The condition of BPS saturation Q3.24j ) implies that 

(sol|if - Z\ sol) = (3.25) 

is valid to any order in perturbation theory. In a sense, H — Z can be viewed as 
a Hamiltonian for the problem with the boundary effects taken into consideration. 
Then Eq. (|3.25|) expresses the complete cancellation of the quantum corrections between 



bosons and fermions. It relates the soliton mass to the central charge, 

M = (so\\H\ sol) = (so\\Z\ sol) . (3.26) 



2 Actually, the supermultiplet is reducible: the states ( | solft) =b | sol / ) ) / \/2 present two one-dimensional 
representations of superalgebra. 
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3.4 Ultraviolet regular izat ion by higher derivatives 

We use higher derivatives to regularize the theory in the ultraviolet. The main advan- 
tage of the method is the explicit preservation of supersymmetry at every step. To 
simplify the construction of the canonical formalism we introduce only spatial deriva- 
tives for the regularization, the Lorentz noncovariance of the procedure does not lead 
to a problem. 

In regularized form the action ( |2.5| ) becomes 



S = iJ d 2 ed 2 x^~D a $(l-^D a $ + W(®)} , (3.27) 

where M r is the regulator mass. In the component form we arrive at the following 
Lagrangian 



H' w ( i -f)* , * + *( i -f) iw+F ( i -f)- 



+2W((f))F- W"((f))ipip} . (3.28) 

The auxiliary field F can be eliminated from the Lagrangian by virtue of the following 
formula: ^ 

F = -( 1 -§?) W - (3 ' 29) 
In perturbation theory near the flat vacuum the bosonic and fermionic propagators 



are 

Ml 



(bosonic) 



M?+pl p 2 - (W^) 2 [1 + (p 2 /M r 2 )]- 2 ' 

M 2 1 
— —t— - — - r , (fermionic) (3.30) 

where Wq is the second derivative evaluated at the vacuum value of the field (p. The 
regularized propagators acquire an extra factor M 2 j{M 2 +p 2 z ) making the integrals over 
p convergent. Strictly speaking, for the quadratically divergent integrals a higher power 
of the same factor is needed - one can achieve this simply by increasing the power of 
(1 - d 2 z /M 2 ) in Eq. Q). 

Proceeding with canonical quantization in the usual way we get the canonical com- 
mutation relations: 



<f>(t, z), (l - ^ J fa z>) 



id {z - z') 



Mt, z), z')\ = ( 7 °) Jb-zf), (3.3i; 
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and the regularized Hamiltonian density Tl, 



n 



1 

2 



ill 2 

r 



1 

+ 2 



3 2 \ 



M 2 



+ d z 



1 - 



M? 



4>-C 



d z (j) - WF - iip 1 - 



M 2 



(3.32) 

7 v+ wvvl • 



It is also simple to construct the regularized conserved supercurrent, 



1 - 



M 2 



flip) + iFd z ftip-i (d z F) fop} 



(3.33) 

Using the canonical commutation relations presented above one can calculate the com- 
mutator of this regularized current with the regularized supercharge. This is a relatively 
straightforward although tedious exercise. In this way we get regularized expressions 
for the energy- momentum tensor and the topological current We will analyze 
them to get the anomaly in the next subsection. 



3.5 Calculation of the anomaly by higher-derivative regular- 
ization 

The particular anticommutator which determines the regularized topological current £ M 
is 



C = -~Tr({^,0} 



7' 



(3.34) 



where J M is presented in Eq. Q3.33|) . The complete expression is rather cumbersome, 
therefore we will limit ourselves to the time component (° which is sufficient for ex- 
tracting the anomaly, 



C° = dAW + 



2M3 



(d z( p)d z F-(d 2 z <P)F 



(3.35) 



where F is given in Eq. ( |3.29| ). The term d z W on the right-hand side is the classical 
topological charge density. The additional terms are due to the regularization, for- 
mally they vanish as 1/M 2 in the limit of the large regulator mass. However, the loop 
integration yields M 2 in the numerator, so a finite additional term survives. 

Using the background field decomposition ( |3.12| ), — > O + X-> we expand the 
right-hand side of Eq. ( |3.35| ) up to the second order in the quantum field x an d keep 
only those terms which survive in the limit M r — > oo, 



2M 2 



(d z <P)d z F-(d 2 z <P)F 



— W 

M 2 



d zX 1 




(3.36) 



Strictly speaking, the Green's function of the quantum field x depends on the back- 
ground 0o; see Eq. (|3.30|) . This dependence drops out, however, and the result for the x 



17 



loop does not depend on the choice of the background field o since the loop integral is 
saturated in the domain of virtual momenta p z ~ M r . Other domains are irrelevant in 
the limit M r — > oo. Thus we can use the bosonic propagator from Eq. fl3.30| ), neglecting 
the (Wq) 2 term in its denominator, 

(^x(l-§) *")=i/(^rf^ [1 + ( ^ /M?)]J =§■ (3-37) 



As a result, combining Eqs. Q3.35|) , ( |3.36|) , ( |3.37|) , we arrive at 



C° = 5 2 |w+^|, (3.38) 

which exactly coincides with Eq. (|3.16|) in the time component. Dimensional regulariza- 
tion and that by higher derivatives lead to one and the same expression for the anomaly 
in the central charge. 

Instead of computing the commutator ( |3.34| ) we can suggest an alternative procedure 
leading directly to Eq. (|3.35| ). To this end it is enough to note that the regularized 
Hamiltonian density given by Eq. ( |3.32| ) can be identically rewritten as follows 

w = C° + JWi-|)^ + (^ + f)(i-|)(W + f) 

j 1 ip + W"ipA , (3.39) 




where the topological charge density (° is the same as in Eq. fl3.35| ) and F is defined 
by Eq. O. 

In the subsequent sections, the very same anomaly will be confirmed, additionally, 
by two alternative analyses: of the effective one-loop superpotential; and models with 
an extra superfield. The extra superfield introduced in the framework of softly broken 
extended supersymmetry serves as a Lorentz-invariant regulator when its mass is taken 
to be large. 

3.6 One-loop nature of the anomaly 

The anomalous term W'/Att was obtained at one loop. A natural question is what 
happens at higher loop orders. In this short subsection we will argue that the expres- 
sion ( |3.16| ) is exact - no higher loop corrections emerge. This should be understood as 



an operator relation. 

The anomaly comes from the ultraviolet range, as we have demonstrated by explicit 
calculations above. In this range one can consider the superpotential term in the ac- 
tion ( ^.5| ) as a perturbation to the first (kinetic) term, which describes massless free 
fields. The superpotential W is proportional to the mass m, and m/M r is our small 
parameter. (In the polynomial model the coupling constant A ~ mv .) 
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Consider, for instance, the operator 7 M J M . At zeroth order in W this operator 
vanishes, in the first order in W it is proportional to W (more precisely, W and its 
derivatives). All higher loops, i.e. higher orders in W, will result in an extra power of 
m/M r as far as contributions from ultraviolet domain are concerned. Thus, the first 
loop is the only one which produces the anomaly. 

Does this mean that if one calculates the Feynman graphs for 7 M J M the only graph 
that contributes is one loop? The answer is no. Certainly, multi-loop contributions 
are present too, but they all represent the matrix element of the operator on the right- 



hand side of Eq. ( 3.18 ). The fact that all nonvanishing multi-loop corrections are of 
infrared origin and, thus, must be referred to the matrix element of the appropriate 
operator, is obvious on dimensional grounds. Thus, the situation with the anomaly 
in two-dimensional models under discussion is a perfect parallel to that in SUSY glu- 
odynamics in four dimensions where the three "geometric" anomalies, considered as 
operator relations, are one- loop exact; all higher-order corrections come as matrix ele- 
ments of the operator TrP^ 2 in the given background field p6| . 

Confirmations of the one-loop nature of the anomaly are provided by calculations 
in the SSG model presented in Sec. |7|. 

3.7 Infrared regularization and fermion-boson cancellations in 
local quantities 

To discretize the spectrum of modes let us put the system in a large box, i.e. impose 
boundary conditions at z = ±L/2. It is convenient to choose them in a form which is 
compatible with the residual supersymmetry ( p.20| ). The boundary conditions we are 
imposing are 

(d z <P + F)\ z=±L/2 = 0, 
(d z + F')^ 2 \ z=±L/2 = 0, 

lkL±L/2 = 0. (3.40) 



It is easy to verify the invariance of these conditions under transformations ( |2.20 ). 



They are also consistent with the classical solutions, both for the flat vacuum and for 
the kink. In particular, the soliton solution O satisfies to d z <p + F = everywhere, and 
boundary conditions (|3.40| ) do not deform it; we can place the kink in the box. 

In Sec. |3.3| the exact equality between the soliton mass and the expectation value of 
the central charge was discussed. Here we will consider a stronger statement, namely 
the exact equality between the energy density inside the soliton and the expectation 
value of the central charge density. This equality is valid to any order in perturbation 
theory and holds locally, for any z. 

Let us use the expression ( |3.39|) for the energy density 7i in the following form, 



n = C° + \ {(0) 2 + (d z <P + Ff + # 2 (<9 2 - + iA(d z + F')4j 2 } , (3.41) 
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where ^1,2 axe the components of the fermion field ip a . We have omitted for simplic- 
ity the 1/M 2 terms on the right-hand side as they do not change the consideration. 
Equation ( |3.41| ) is the exact operator form. In the classical limit, where ip = 0, the 
expression is positive-definite. Its minimum is achieved provided the field <fi is chosen to 
be time-independent and to satisfy the BPS condition d z cf) + F = 0. Then the soliton 
mass is given by the integral over the first term in Eq. (|3.41| ), i.e. by the central charge 
Z. But now we want to analyze the relation (|3.41| ) in its local form. 

For the classical kink solution the difference 7i — C° obviously vanishes at any z. 
This vanishing persists at the quantum level, 

(sol\H(x) -C°(aO|sol) = 0. (3.42) 

Indeed, the operator Tt(x) — (°(x) can be presented as the anticommutator |T^ 

{Q 2 ,J°(x)} = H(x)-(°(x), (3.43) 

which vanishes upon averaging, since Q%\ sol ) = 0. As a result, not only do the quantum 
corrections to the kink mass coincide with those for the central charge, but the energy 
density distribution coincides with the central charge density to all orders. 

It is instructive to check in more detail how the cancellation works in the one-loop 
approximation. To this end we expand the bosonic field near the classical soliton 
background (fi Q , 

<P = <p + X- (3-44) 

In the quadratic in the quantum fields x approximation th difference TC — £° takes the 
following form: 

H - C°] quad = ~ {x 2 + [(d z + F') X f + iih(9z - F')iPi + iMd z + F')i, 2 ) , (3.45) 

where F' is evaluated at <fi = <p and the subscript "quad" denotes the quadratic ap- 
proximation. We recall that the prime denotes differentiation over </>, 

dF _ d 2 W 

= d(j)~~ d(f) 2 ' 

From this expression one can see that a convenient basis for the mode expansion of x 
and ip2 is given by the eigenfunctions Xn{z) of the Hermitian differential operator L 2 , 

L 2 xM=<4Xn(z), L 2 = P^P, (3.46) 

where the operators P and P* defined as 

P = d z + F' , P^ = -d z + F' (3.47) 

are Hermitian conjugates. The boundary conditions for the modes in the box follow 
from Eq. ( gg ), 

(d z + F')xn(z = ±L/2)=0. (3.48) 
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As for the basis of the mode expansion for ipi, it is formed by the eigenmodes of the 
operator L 2 , 

L 2 x n (z)=u 2 n Xn(z), U = PP\ (3.49) 
with the boundary conditions 

Xn{z = ±L/2) = . (3.50) 

With boundary conditions (|3.48|) and (|3.50|) all eigenvalues of the operators (|3.46|) 
and ( |3.49| ) are the same, with the exception of the zero mode xo ^ d(j)o(z)/dz in 
Eq. ( |3.46|) . The operator ( |3.49| ) has no zero mode. Moreover, the eigenfunctions Xn 
and Xn are algebraically related, 

Xn = —PXn, Xn = —P ] Xn- (3.51) 

The expansion in eigenmodes has the form, 

X(x) = J2 b n(t) Xn(z) , 1p 2 (x) = J2 Vn(t) Xn(z) , 

^i{x) = Y J Ut)Xn{z). (3.52) 

Note that the summation does not include the zero mode Xo{z)- This mode is not 
present in ipi at all. As for the expansions of x an d ip2, the inclusion of the zero mode 
would correspond to a shift in the collective coordinates z and r\ (see Sec. |2~2| ). Since 
we deal with the bosonic kink centered at z — 0, these shifts must be forbidden. 

The coefficients a n , r\ n and £ n are time-dependent operators. Their equal time 
commutation relations are determined by the canonical commutators (|3.3|), 



\Pmi bn\ i$mn > \jlmi Vn\ $mn > {'Cm; ^n} &mn • (3.53) 

Thus, the mode decomposition reduces dynamics of the system under consideration to 
quantum mechanics of an infinite set of supersymmetric harmonic oscillators (in higher 
orders the oscillators become anharmonic). The ground state of the quantum soliton 
corresponds to setting each oscillator in the set to the ground state. 

Constructing the creation and annihilation operators in the standard way we find 
the following nonvanishing expectations values of the bilinears built from the operators 
a n , r\ n and £ n in the ground state: 

(^n)sol = Y , (^)sol = , (^n^n)sol = ~ • (3.54) 

The expectation values of other bilinears obviously vanish. Combining Eqs. Q3.45| ), 
( jT5g ) and ( |334D we get 
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£ { y x1 + 2^ [{dz + F ' )Xn]2 - y x1 ~ ir [{dz + F ' )x A s ' (3 - 55) 



The four terms in the braces in Eq. (|3.55| ) are in one-to-one correspondence with the 



four terms in Eq. (3.45). Note that in proving the vanishing of the right-hand side we 



did not perform integrations by parts. The vanishing of the right-hand side of (|3.45 ) 



demonstrates explicitly the residual supersymmetry (i.e. the conservation of Q 2 ) at 
work. The cancellation occurs in the very same manner with a finite regulator mass 
M r due to the residual supersymmetry. 

Equation (|3.42|) must be considered as a local version of BPS saturation (i.e. con- 
servation of a residual supersymmetry). The nonrenormalization of the expectation 
value of 7i — C° over the soliton is a direct analog of the nonrenormalization theorem 
in the instanton background leading to the exact (3 functions |16|] in supersymmetric 



gluodynamics. The idea of fermion-boson cancellation of the quantum corrections in 



local quantities was first formulated in Ref. |10|. However, in this work the problem 
of ultraviolet regularization was ignored. Consequently, the quantum anomaly in the 
central charge density was lost. 

In general, the phenomenon of fermion-boson cancellation takes place for the expec- 
tation value of every local operator which can be represented as {Q 2 , 0{x)} where 0(x) 
is some fermionic operator. We will use this feature later on to calculate the quantum 
corrections to the energy distribution and to the mean field in the kink. 

3.8 Explicit expressions for the modes and Green's functions 

In this subsection we present explicit expressions for the modes and Green's functions 
in the SSG model. A nice feature is that everything is expressible in terms of the modes 
and Green's functions in the flat vacuum. A similar relation exists for the SPM model, 
and it is outlined at the very end of the subsection. Green's functions in the soliton 
background are needed for the calculation of the field and energy profiles at one loop; 
see Sec. |5|. 

Let us consider the SSG model. In the kink background the function F'(z) which 
determines the operators P and P^ defined in Eq. (|3.47|) has the form 



F' = m tanh mz , (3.56) 
Consequently, the operators L 2 and L 2 which determine the mode decomposition are 

L 2 = pt P = + m 2 _ 2 r f (3.57) 

cosh mz 

and 

L 2 = PP ] = -dl+m 2 . (3.58) 

The key observation is that the operator L 2 is a quadratic operator in the flat vacuum, 
and its eigenmodes are those of free motion, 

Xn = smp n (z + — ^ . (3.59) 
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These solutions satisfy the boundary conditions ( |3.50|) at 

p n = — : n = l,2,3,..., (3.60) 

with the eigenvalues 

u 2 n = m 2 +p 2 n . (3.61) 

As discussed in the previous subsection, see Eq. ( |3.51|) , all nonzero modes Xn of the 
operator L 2 are algebraically expressible via modes Xn- The two operators, L 2 and L 2 , 
form a superpair in the sense of Witten's supersymmetric quantum mechanics [2^,|27] 
(for a pedagogical discussion see p8|). All non-zero modes of the operator ( |3.57|) are 
then trivially calculable by virtue of the algebraic relation ( p. 51 ). 



Xn{z) 



LO r , 



[—d z + m tanh mz] sin 



Pn[Z + 



L 



A similar relation holds for the Green's function, 



Xn{zi)Xn{z2) 



0J' 



(3.62) 



(3.63) 



G u (Zl, Z 2 ) 



-d zi + m tanh mz\){d Z2 + m tanh mz 2 ) G^(zi, z 2 ) — ^=0(^1; z 2j 



(3.64) 

where G^zi, z 2 ) is the Green's function for the free operator L 2 . In the limit of a large 
box, L 00, the relevant expression for this Green's function is 



G LU (z 1 ,z 2 ) = _ ^: ex P (~ k ~ z ?\) 



(3.65) 



with k = \/m 2 — uo 2 . Note that Eq. ( p.64|) produces the Green's function orthogonal to 
the zero mode of the operator L 2 . We will make further use these relations in Sec. [5| 
and H 



In the general case, Witten's superpair of operators has the form (see e.g. [28 



o2 2 K{K-\)m 2 

d 2 z +m 2 and 

cosh mz 



d 2 + m 2 



K(K + l)m 2 
cosh 2 mz 



(3.66) 



where K is integer. We used this fact above at K = 1, to relate the eigenmodes in the 
SSG kink background to those of free motion. We can now repeat the trick to relate 
the eigenmodes in the SPM kink background to those in the SSG kink background, 
explicitly written above. Indeed, for the SPM kink, 



F'(z) = m tanh 



mz 



(3.67) 



and the quadratic operator which determines the mode decomposition for the field x is 

1 



ptp 



-dl + m 2 



2 cosh(mz/2) 



(SPM) . 



(3.68) 
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The sister operator is 



i-i 



2 cosh(mz/2) 



(SPM) . (3.69) 



It is not difficult to see that 



Lf^z) = \Lf G {z/2) + \m\ (3.70) 
3.9 Separating soliton from the boundaries 

We pause here to emphasize one important point concerning calculations with bound- 
ary conditions. The relation between the eigenmodes in the two sectors connected by 
the quantum-mechanical supersymmetry operator P and the exact coincidence of the 
spectra (with one extra zero mode in the kink sector) only holds when the boundary 
conditions in the two sectors match each other. Namely, these conditions at z = ±L/2 
should be the same for Pxn/^n in the kink sector as for Xn in the dual sector (the 
latter is the free motion sector in the case of the SG model[|). Under these conditions 
we found the exact local in z boson-fermion cancellation of the difference between the 
energy density and that of the topological charge in Eq. (|3.55 ). However in other quan- 



tities, e.g. in the density of the topological charge, to be considered in the next section, 
as well as in non-supersymmetric models (like the SG model discussed in Sect. |6|), the 
local cancellation does not take place. In calculation of such quantities one should 
clearly realize the need for separating the effects of the soliton background from those 
introduced by the boundaries. The criterion for the separation is obvious: the effects 
of the boundaries are localized near the boundaries (at distances set by the mass m 
and by the regulator parameter M r ) and are carried away to infinity in z in the limit 
of infinitely large box: L — > oo, while the effects of the soliton proper stay at finite 
distances, whose scale is set by the mass m. This behavior is illustrated in Fig. [I], 
where the quadratic average of the quantum fluctuations, x 2 { z )) is shown for the SG 
kink background and for the flat background, as calculated with discrete modes in a 
finite box with a finite ultraviolet regulator parameter. The spikes near the edges of 
the box are the discussed artifact of the boundary conditions and are moved to infinity 
in the limit of an infinite box. The soliton background effect is the difference between 
the curves at the center of the plot, and it remains fixed in the limit L —>■ oo. 

Thus the effects of the soliton proper can be separated by considering the system 
in a large box up to distances well inside the box, i.e. up to z = ±(L/2 — A), where 
it is sufficient to choose A such that L ^> A ^> mT 1 . This interval of z where the 
soliton effects already have approached their asymptotic behavior, but the effects of 
the boundaries of the box are still not essential, can be refered to as fiducial. In what 



3 This is not the procedure used in the standard calculations so far: rather the same (usually periodic) 
boundary conditions are imposed on \ n and \ n (see e.g. in p9[), which introduces a difference in the 
spectra through a phase shift. 
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Figure 1: The quadratic average x 2 { z ) i n a finite box with a finite ultraviolet regulator 
parameter in the SG kink background (solid line) and for a free field (dots, coalescing 
into thick line away from the edges of the box). The arrows show the boundaries of the 
fiducial interval, that encloses the effects of the soliton proper as opposed to those of 
the boundary conditions. 



follows we will imply that only the quantities inside the fiducial interval are considered 
and the integral quantities are also evaluated over the fiducial interval. 

In practice the extraction of the quantities related to the averages over the quantum 
fluctuations within the fiducial interval is readily done by choosing the appropriate 
Green's functions for the fields of the fluctuations. Namely, in the limit of the box and 
the fiducial interval both taken to infinity these Green's functions should be chosen with 
the Feynman boundary conditions, i.e. they vanish when one of the arguments z goes to 
either plus or minus infinity, provided that u is shifted to the upper complex half-plane. 
This is the boundary condition under which we have specified the free-motion Green's 
function in Eq. ( 3.65| ). The discussed above quantum-mechanical supersymmetry then 
generates the appropriate Green's functions in the soliton sector of the SG model as 
well as of the SPM model (cf. Eqs. (|3.64|) and ( |3.70| )), which Green's functions we will 
be using in the next section. 
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4 The one-loop correction to the kink mass 



4.1 Central charge and soliton mass at one loop 

Once BPS saturation is proved the actual calculation of quantum corrections to the 
kink mass becomes a simple job. Indeed, according to Eqs. ( |3.26|) and (|3.17|) 

M = (sol|Z|sol) = (sol|W + W"|sol> 2=+00 - (sol| W + W"|sol),=-« , (4.1) 

47r Air 

where the term W"/47r is due to the anomaly we discussed above in detail. At the one- 
loop level this term can be substituted by its classical value at z = ±oo. Additionally, 
all we need to know is the one-loop quantum corrections to (sol | W| sol) at z — > ±oo. 
These corrections are the same as for the flat vacuum where the standard perturbative 
expansion can be used. 

At z — > +oo the classical kink field 4>o goes to the constant 0+. In the flat vacuum 
the operators and W can be presented as 

4>{x) = <pt + xO) , w[4>{x)) = w + \w£ X 2 (x) + o( x 3 ) , (4-2) 

where Wo = W[0+] and W ' = W" [</>+] are c numbers and x( x ) is the operator. In the 
one-loop approximation 

(W) = Wo + \ W ' ( X 2 (x)) q = Wo + \ W ' G(x, x) , (4.3) 

where G(x, x') is the propagator of the field x- Using the regularized momentum trans- 
form of Eq. ( |3.30|) we get 

2 v r d 2 p M 2 r i = J_ 4M r 2 

X /0 ~ J (27T) 2 M r 2 + P 2 p 2 _ [! + (p2 /M 2)]-2 - 4?r 11 (W ')2 ' 1 ' j 

where M r is the ultraviolet cut off. Equations ( |4.3| ) and ( |4.4| ) give the one-loop result 
for (W(z — ► +00)). The result for (W(z — > — oo)) is different only by the sign. 
Thus, for the kink mass we get 



W n ' 

47T 



1 4M r o' 



(4.5) 



The first term is the classical result for the mass, the one-loop correction is given by 
the second term, the logarithm in the square brackets arises from the correction to (W) 
and the additional term 2 is from the anomaly. 

It is clear that to make sense of the quantum correction we need to formulate the 
renormalization procedure, i.e. to fix some observables other than the kink mass. The 
mass of the light particles proves to be appropriate. It is simple to evaluate the one-loop 
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correction to the fermion mass in the SSG model. It is given by the tadpole diagram 
and can be written as 

m ren = - (W"> = -Wq ~ \ Wr ( X 2 (x)) o , (4.6) 

where, as above, (x 2 ) appears. At the one-loop level in the SSG model m ren is both, 
the on-shell mass, and the renormalized parameter in the superpotential (see the next 
subsection). Accounting for the fact that in the SSG model W = mv 2 smcj)/v and 
4"t = the result (4.5) for the kink mass in terms of m ren has the form: 



M SS G = 2m r 



2 1 

4:71 



2m rcn v 2 -^p. (4.7) 

Z7T 



The logarithm term is absorb by the mass renormalization, what remains is the anomaly 
term. 

In the SPM model the one-loop expression for the on-shell mass is more complicated, 
it will be obtained in Sec. f).2\ We can introduce, however, an off-shell renormalized 
mass parameter as follows: 

Mren = ~ (W) Q = 2A vac = 171 + 2A ( X )n . (4.8) 



Using the fact that (W) = 0, see Eq. (|2.27|) , we can express (x) in terms of (x 2 )o 
which was determined above, 



Then, 



M S pm = ^-^. (4.10) 



1 /^ren _ Mr 

6 A 2 2tt 

The definitions ( |4.6| ) and (|4.8|) of the renormalized mass parameter are the same as 
those adopted in Refs. [p^,|I^|, and our results coincide. What then is the difference 
between our calculations and those of our predecessors? 

What was correctly calculated in Refs. fl4], 15], as well as in some previous works, 



was (sol | H | sol). The problem was that the result did not match (sol|^|sol). What we 
present here is the calculation of (sol 1 2 1 sol). We resolve the problem by adding the 
anomalous term in Z. 

The next section presents an effective Lagrangian approach to the calculation of 
Z which reproduces the very same anomaly. Moreover, in SPM we suggest another 
definition of the renormalized (off-shell) mass parameter, which is most natural from 
the point of view of the effective Lagrangian. 



4.2 Soliton mass from effective superpotential 

A convenient notion to use is the effective superpotential W e ff(0) which appears in 
the background field approach to perturbation theory. It accounts for the quantum 
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corrections in the potential energy which is considered as a function of the background 
field bck . 

Let us calculate one-loop corrections to W e ff(0) choosing as the background constant 
fields and F, 

0(x) = 0bck + x(x), F(x) = F hck + f(x). (4.11) 

The difference with Eq. (|4.2| ) is that the background field is not a solution of any classical 
problem and can be chosen arbitrarily. Then the part of the Lagrangian ( |2.7| ) quadratic 
in the quantum fields is 

£( 2 ) = - [d^X + f + 2W"fx + FW"'x 2 + $i ~ WW) . (4.12) 
The one-loop correction to the effective Lagrangian £ e ff(0bck> -^bck) is defined as 

exp^y d 2 xA£ cfr } = JvxVtfjexp^iJ d 2 a;£ (2) j . (4.13) 



As is seen from Eq. ( p.7|) it is sufficient to determine the part of A£ e ff linear in in 
the limit Fb C k — > to get the one- loop correction to W' cS , 



AW' cS 



W" , 4M 2 
In ■ 



■ill 



(W") 



(4.14) 



- ^bck 



As mentioned above, the difference with the previous calculation of (x 2 )o is that 0bck 
is a parameter not equal to its classical value 4>t • 

It is easy to construct the one-loop correction to the superpotential WW from its 
derivative (|4.14|) . The total WW is 



W, 



cff 



W 



8vr 



4M 2 
ln^Z^ + 2 



n\2 



(4.15) 



Note the appearance of the additional 2. This is another alternative way to reproduce 
the anomaly in the central charge. The origin of this 2 can be traced back to the 
differentiation of ln0 in the superpotential. Thus, within the effective Lagrangian 
approach the anomaly emerges as an infrared, rather than ultraviolet, effect. This is a 
usual story - both faces of the anomaly, infrared and ultraviolet, are intertwined. 
The condition W' eS = defines the shift of the mean vacuum field, 



0vac = (vac|0|vac) = 0+ + 01 , 0] 



1 Wg , 4M3 
tt m 



8tt Wg (W£) 2 



(4.16) 



At this value of the field 

W eff (0vac) 



87T 



4M 2 

ln (w# + 2 . 



(4.17) 
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The kink mass is denned by the difference between W e g at z = ±00, 

M = 2W eff (0 vac ), 



(4.18) 



which coincides with the result ( |4.5| ). Thus, we rederive the anomaly from the effective 
Lagrangian. 

Let us stress that the approach presented above is based heavily on supersymme- 
try. Although we did not use the superfield formalism explicitly it was implicit when 
we introduced the F component of the background field independent of (p. Had we 
eliminated the F term from the very beginning the effective Lagrangian would differ by 
finite (nonlogarithmic) terms. 

The effective Lagrangian approach is also convenient for renormalization, all ul- 
traviolet logarithms can be hidden in the parameters of W e ff. In particular, in the 
polynomial model 



m 



W cff = — < 



A 



A 

47/ 



M 2 
ln_i_ + 2 



A 2 2 



4A r 3 

It is natural to define the renormalized mass m ren as 

A 2 



m 

2A Airni 



A , 4M 2 
In — ^~ 



(4.19) 



m, 



effl 



m — 



2-nm 



(4.20) 



In terms of m ren the expressions for W e g and 



m. 



A 



off 



4A " 3 
and the soliton mass is 



^ A 1 m ren 

4^ 01n 4W' 



contain no ultraviolet cut off, 
A 



m, 



2A 



+ 



27rm, 



(4.21) 



1 TJV^ 

M SPM = - 



(4.22) 



This expression superficially looks the same as the classical expression. Of course, one 
can do things differently, expressing M in terms of /i ren = 2A0 vac introduced in Eq. 
( |L8p . Then we recover Eq. ( 4.10| ). The result for the soliton mass expressed in terms 
of the on-shell mass m of the elementary quanta is given in Eq. ( |5.26| ). 
In the sine-Gordon model the one-loop effective potential reads as 



cff 



mv 2 sm((p/v) — 



m sin((f)/v) 



8tt 



In 



4M 2 



m? sin (4>/v] 



+ 2 



so that the vacuum field gets no shift 0, 



0vac = 4>Q = TTv/2 . 



(4.23) 



(4.24) 



4 In fact, vac 



stays at irv/2 to all orders. 
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while the renormalization of the mass parameter is given by 



m ren = m - In — — . (4.25) 



1 7r v m' 



This is identical to Eq. ( [4.6[ ). Note that in the sine-Gordon model at one loop m ren in 
Eq. ([4.25 ) is the same as the on-shell mass fh of the elementary quanta (see Sec. §). 



The soliton mass is the same as in Eq. ( [4.71 ). 



5 Field and energy profile at one loop 

5.1 General relations for profiles 

The quantum corrections to the kink mass, see Eqs. ( |4.5[ ) , (4.18), are defined by the 
simple dynamics of the flat vacuum at z — * ±oo. Here we will find the one-loop 
corrections to two local quantities: the mean field and the energy density in the kink. We 
start with a general superpotential W(<&), and then demonstrate the explicit solution for 
the polynomial and sine-Gordon models. Similar corrections to the domain wall profile 
in four dimensions were found previously in non-supersymmetric models in [^1|,^2[ and 
in super symmetric theories in ]18| . 
The mean field (f> (z) is defined as 

0(z) = (sol I <f>(t, z) I sol) = <j) {z) +<j) 1 {z) , (5.1) 

where we denote as <pi{z) the deviation of the mean field from the classical kink solution 
4>o(z). To derive the equation for the mean field let us start from the anticommutator 

{A,Q2} = -(d z <P + F) , (5.2) 

where F denotes — W'(0). As was discussed in Sec. |3]7| the expectation value of the 
anticommutator vanishes due to BPS saturation, Q 2 |sol) = 0. Therefore, we get for 
the mean field 

_ = -(soil F\ sol) . (5.3) 
az 

Expanding the operator field <fi near the classical solution, 

<f> = (j> {z) + <t> 1 {z)+x{t, z) , (5.4) 

we get from Eq. ( |5.3| ) 

■^M*) = -F'M') ~ \F"(x\z, t)) sol ■ (5-5) 

Throughout this section the derivatives of the superpotential with respect to the field <fi, 
denoted by the appropriate number of primes, are evaluated at the classical profile, i.e. 
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at (j> — 4>o( z )- One should also notice that the quadratic average of the quantum fluc- 
tuations (x 2 (t, z)) is in fact time- independent, since it is calculated in the background 
of a static field configuration. 

Let us first consider the behavior of the solution of Eq. ( |5.5|) at the infinities. For 
the limiting behavior at z — > ±00 we found in Sec. 01, see Eq. ( J4.16 ), 



lim 0i(2;) 



w, 



2W, 



(X 2 (z, t)) 



*-±oo 8vr Wo (W, 



z— >±oo 



(5.6) 



These limits are consistent with Eq. (|5.5|) , its left-hand side is vanishing at z — > ±00. 

Using the expansion ( |3.52| ) of x i n eigenmodes and Eq. ( ft. 54 ) we find for the 
quadratic average 

,2/ 



X^))sol=E^ 



The same result can also be presented as 

(X 2 (%) )soi = G(x;x) 



(5.7) 



(5. 



where G(x; x) is the limit at coinciding space-time points of the Green's function 
G(x; x'), satisfying the second order equation 

'-d 2 t +(d z -F')(d z + F')] G(x;x')=i5(t-t')[5(z-z')- X o(z)xo(z')} • (5-9) 

The subtraction of Xo( z )Xo{ z> ) m the right-hand side is needed for the exclusion of the 
zero mode. 

Due to the time independence of the background, the Green's function is conve- 
niently expanded in frequency modes: 



-00 Z7T 



G{t;z;t',z') = 

where G u (z; z') satisfies the one-dimensional equation 

+ (d z - F') (d z + F')] G u (z; z) = 1 [6{z - z) - Xo(z)xo(z') 



(5.10) 



(5.11) 



The solution to this equation for G^(zi, z^) can be written in terms of the normalized 
eigenfunctions Xn{z) and the corresponding eigenvalues u n defined in Eq. (|3.46|) in the 
standard way, 



G u (z] z 



uj z — ujz + is 



At coinciding points we get 

(X 2 (%) )soi = G(x;x) 



+00 



2^ 



(5.13) 
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and the integration over uj reproduces the representation ( |5.7| ). 

In terms of the quantity (x 2 (z)) the solution to the equation ( |5.5| ) for the correction 
to the field profile can be written as 

Uz) = A Xtt{z) f ^hmm iu = j F(Mz)) f^mM> Au , 

2 Jo xo[u) 2 Jo F(<p (u)) 

(5.14) 

where the last equality makes use of the fact that the translational mode Xo is pro- 
portional to F due to the BPS saturation at the classical level. The lower limit in the 
integral is the collective coordinate for the soliton center which we have chosen to be 
at z — 0. 

Once (x 2 ( x ) ) an d 4>i are found the quantum correction to the energy distribution 
in the kink is also fully determined. The quantity 

W'(z)\ W" 2 , 



~ /sol 27r 



4tt \ F" / Vl F" dz 



(5.15) 



represents the total energy located in the interval between — z and z (implying that z > 
0). In the last expression W and its derivatives over (ft are evaluated at (f> = (f)o(z). In the 
large z limit M(z) goes to the soliton mass M. The energy density is (l/2)dM(z)/dz. 

Thus, the one-loop corrections to the mean field and to the energy density are 
expressed via integrals over the Green's function G w (z; z) defined by Eq. ( |5.11| ). These 
results were obtained for a general superpotential function W. In the special cases of 
the polynomial and sine-Gordon superpotentials (see Eqs. ( |2.8| ) and ( |2.9| )) this Green's 
function is known analytically and we will get explicit expressions for the one-loop 
corrections. 



5.2 Polynomial model 



The partial Green function G^z^z-i) is readily constructed by twice applying the re- 
lation ( |3.64|) to the free Green's function in Eq. ( |3.65| ) (it can also be read off from 
Refs. f2l],|22|). At coinciding points G w {z\ z) is given by 

3 m 2 / 



Gu(z\ z) 



2 y/m 2 — uj 2 



1 

1 ^ 



uj 2 4 cosh (mz/2) 



771 



1 3 m 2 tanh 2 (mz/2) 
3m 2 — Auj 2 cosh 2 '(mz/2) 



(5.16) 



Thus, substituting this expression in Eq. (|5.13|) and performing Euclidean rotation in 
the integration over u (i.e. the substitution u — > i uj) one finds for the average over the 
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fluctuations 



sol 



+00 



du 



Atv \Jio 2 + m 2 



1 

1 + ^7 



3 m 



u; 2 4 cosh (m;z/2) 



7?). 



4tu 2 + 3m 2 

1 , 4M 2 
— In — — 

47T 



3 m 2 tanh 2 (777,2/2) 
cosh 2 (mz/2) 
3 1 



(5.17) 



+ 



1 tanh 2 (m^/2) 



8ir cosh 4 (mz/ 2) 4 cosh 2 (mz/2) 



It is understood here that the divergent part of the integral is regularized by the ultra- 
violet regularization procedure described above. Upon substitution in Eq. ( |5.14j ) this 
expression gives the final result for the correction to the profile of the field across the 
soliton: 




. 4M 2 2tt\ 
In — f + -= 
m 2 y/3j 

tanh(mz/2) 

4 TV J cosh 2 (mz/2) 



d<fic 



dm 



A 



(5.18) 



The first two terms in this expression clearly can be absorbed in a redefinition of the 
parameters m and A in the classical profile of the field. Thus the full one-loop corrected 
field d> across the soliton can be written as 



4> = M z ) + M z ) 



m A 

— tanh(mz/2) + 

2A 771 



1 



3 \ tanh(mz/2) 
i 2\/3~ r 47r/ cosh 2 (m2/2) 



where 



A 5 



777 



27r m 



In — + — ) and A = A 



A a 



777 z 



(5.19) 



(5.20) 



are the appropriately renormalized parameters. It should be noted that the mass pa- 
rameter in this renormalization scheme exactly corresponds to the on-shell mass renor- 
malization. In other words, the physical (pole) renormalized mass of the scalar boson 
(and, by supersymmetry, of the fermion) is equal to m, as can be seen from the expo- 
nential approach of the profile in Eq. fl5.19|) to the vacuum value at z — > 00. 

Substituting the result fl5.18Q for the profile correction into Eq. (|5.15|) we find the 
energy M(z) sitting between — z and z, 



M(z) 



9 

m 



6A Z 



tanh(777z/2) 



m 




tanh (7772/2) + 777 

2 TV 

The first line is similar to the classical energy distribution up to the transition to 
renormalized parameters. The correction factor in this line is equivalent to an additional 
finite renormalization of A. 
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The qualitative feature which quantum corrections bring in is that M(z) approaches 
its asymptotic value at large z more slowly. Classically (the first line) it approaches as 
exp(— 2mz). The one-loop correction contains exp(— mz). This can be simply under- 
stood in perturbation theory: the energy density operator can only produce one-particle 
states at the quantum level. 

We also find it quite instructive to consider a "kink effective superpotential" W e fr(0), 
which generates the full one-loop profile of the field, as the solution to an "effective" 
equation 

^0 = W eff (0) . (5.22) 

This effective function can readily be found by noticing that the expression for the 
average (x 2 (- 2; ))soi given by Eq. ( |5.17| ), which enters the quantum BPS equation ( |5.5| ), 
is in fact a fourth power polynomial in the classical solution <f)o(z). Thus the inhomo- 
geneous term in Eq. ( |5.5|) can be written as the derivative of a fifth power polynomial 
in 0. Hence one arrives at the following expression for W c g 



I7T 



A 2 —3 ^ 16A 4 —5 



3 m 



5 m 



(5.23) 



It is important to emphasize that the so defined soliton effective (super)potential 
WW is different from the Coleman- Weinberg type effective (super)potential W c s in 
Eq. ( |4.15|) . The latter describes the response of the action to a space-time constant 
field, and is appropriate for finding the quantum corrections to the vacuum field, while 
the purpose of the former is to describe the quantum corrections to the soliton profile 
through the equation ([5.22|) . In fact, the W e ff also describes the correction to the vac- 
uum field as the asymptotic value of the solution at z — > oo, which perfectly agrees with 
vac found from W c s: W cff (0 V ac) = 0. In addition the soliton effective superpotential 
encodes all information about the detailed behavior of the field profile at finite z. In 
particular, its second derivative gives the on-shell mass of the elementary particles, 



W, 



cff 



-m . 



(5.24) 



unlike the mass parameter m ren , calculated in Eq. ( [4.20|) from the curvature of W c s 

X 2 



m = m ren + 



1 



7T 



Tcm 



V3 ■ 



(5.25) 



Thus, the mass of the kink from Eq. (|4.22|) , being expressed through the on-shell mass 
parameter m reads as 



1 m 3 m m 1 m 3 m m 
= 6A^ + 273~2T = 6X 2r+ 6V3~2T 



(5.26) 



The same expression for M appears in Eq. ( |5.22| ) as the asymptotic value of M(z) at 
large z. The first two terms in Q5.26D are equivalent to the result of Ref. |L0| where it 
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was obtained within a similar approach. However our result differs by the third term, 
which arises from the anomalous part of the central charge. 

A remarkable property of the soliton effective (super)potential is that it is a finite 
polynomial of the field, unlike that of the Coleman- Weinberg type. This is due to 
the phenomenon of "nullification" [ZI,|22], 3"0"|: the vanishing of all on-shell amplitudes 
for multiboson production at threshold, starting from a finite number of the produced 
bosons. The on-shell amplitudes are related to the profile of the soliton field |21], RT| . 



5.3 Sine-Gordon model 



Following the same routine as for the polynomial model, we write the expression for 
the partial Green's function G UJ (z; z) in the background of the classical soliton profile 
given by Eq. ( [2.12|) , 



Geo (z; z) 



2^/m? 



10* 



m 



1 



lo 2 cosh 2 mz 



m" 



LO 



rn 



(5.27) 



The average over the quantum fluctuations then takes the form 



sol 



duo 



1 . AM 2 
m — -f- 



oo 4 7r V 'lo 2 + m 2 
1 1 



i + 4 



m 



m* 



lo 2 cosh mz 



m 



... , , . (5.28) 

An m 2 2tc cosh mz 

Thus, according to Eq. Q5.14|) , the one- loop correction to the soliton profile is given by 
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(5.29) 



Absorbing the first term in this expression into a renormalization of the mass in the 
classical field profile, one can write the final expression for the full one-loop profile of 
the soliton field <b as 



2v arctan[tanh(mz/2)] + 



1 tanh mz 



Anv cosh mz 



where 



m 



m = m — 



. AM 2 
m — ~~- 



(5.30) 



(5-31) 



5 m* 

is the renormalized on-shell mass of the elementary particles in the sine-Gordon model. 
One can readily notice that in this model at one loop the renormalized pole mass 
coincides with the one defined from the one- loop effective potential: m ren in Eq. ( |4.25[) . 
This is no surprise, since in this model at one loop there is no dispersion of the particles' 
self-energy. (This property, and the relation m = m ren is lost, however, at the two-loop 
level, see Sec.0.) 
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The amount of energy located in the interval between —z and z is readily calculated 
from Eq. ( |5.15| ) and is given by 



, . / 9 1 \ , m tanhmz 

M(z) =2m (v 2 - — tanhmzH - 1 — . (5.32) 

V 4:71/ n cosh mz 

One can notice that in this case the approach of M(z) to its asymptotic value M is 
described by terms proportional to exp(— 2mz), i.e. corresponding to propagation of a 
two-particle state. This is due to the fact that in the sine-Gordon model the number of 
particles is conserved modulo 2, so that the energy operator does not couple through 
quantum loops to a one-particle state. 

The soliton effective superpotential is also readily recovered in this model and reads 
as follows, 

VW — 2 • 4> . m ( ■ 4> 1 • 3 0\ /, q ^ 

VV c s = mv sin — h - — sin sm — I , (5.33) 

v A Ti \ v 3 v J 

and is polynomial in sin (<p/v). 



6 "Supersymmetry" of the non-SUSY sine-Gordon 
model. 

Here we will discuss a "supersymmetry" in the problem of calculating the soliton mass 
(at one loop) in the non-supersymmetric sine-Gordon model, whose Hamiltonian is 
simply the bosonic part of that of the SSG model. The sine-Gordon Hamiltonian 
density can be written as 



n = d z w+ l - {0 2 + (<9 Z + Ff) (6.i) 



with W(0) given by Eq. ( |2.9| ) and F((f>) = — W(0) = — mv cos(0/t>). Since the 
purpose of this section is mainly illustrative, we will ignore the necessity for ultraviolet 



regularization. This can be done along the lines of Sec. |3.5| . We will not introduce 
the explicit regularization to avoid cumbersome formulae which could overshadow our 
key points. One should understand that in doing so we are bound to miss some finite 
(non- logarithmic) piece in the soliton mass. This flaw could easily be eliminated, if 
desired. 

In order to write the one-loop (quadratic) expression for the energy in terms of the 
quantum fluctuations x over t ne kink background <j)o{z) (see Eq. ( f2.12| )) one should 
take into account that, unlike the supersymmetric case, the flat vacuum energy density 
receives a one-loop renormalization and is to be subtracted. The proper expression is 
as follows: 



1 



2 



(^)kink - W vac -d t W=-\x 2 + [d g + F) X ~ X ~ {dzX) - ™ 2 X 2 • (6-2 
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Here F' = m tanhmz is the background value at 4>o{z) ) and W is the profile of 
mv 2 sm((j)/v) across the soliton (including the one- loop correction). Finally, x de- 
notes in this section the field fluctuations over the flat vacuum at O — nv/2 (see Sec. 



We consider the contribution of the quantum fluctuations x an d X, explicitly written 
in Eq. ( |6.2| ), by expanding these fields in the appropriate eigenmodes. Namely, the field 
X is expanded in the eigenmodes Xn( z ) of the operator 



L 2 = -(d z -F') (d z + F') = -d 2 z +m 2 - 



2 m 2 



cosh mz 



(6.3) 



while the corresponding operator for the modes Xn of the field x m the flat vacuum is 
obviously 

L 2 = -d\ + m 2 . (6.4) 
Since the operators L 2 and L 2 are related by a quantum-mechanical supersymmetry 



(see Sec. |3~8|) , we use here the fact that the spectra of both operators coincide (modulo 
the zero mode of L 2 ). We also use the ensuing relations for the eigenmodes. It is seen 
that the soliton sector can be dubbed "bosonic" while the flat vacuum "fermionic". 
These two sectors are related by a quantum-mechanical "supercharge" P = (d z + F ). 

With our boundary conditions, the terms in Eq. ( |6.2| ) can be formally f] written as 
follows: 
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2 UJn^O Un ^ 2 



n ~2 / \ 
Xn( Z ) 



(6.5) 



Here the relation 



PXn{z) = U n Xn{z) 



is used in the second line. Note also that the total derivatives are explicitly kept. 
Using these relations, the expression ( |6.2|) for the kink energy density can be formally 
rearranged as 



(H) kmk - (H) vac - d z W = --dl £ 



2 Wn^O ^ 



XnO) - Xn( Z ) 



(6.6) 



5 As was mentioned, for illustrative purposes we ignore here the problem of the specific ultraviolet 
regularization of the one-loop terms in the sine-Gordon model. Thus, an analog of the anomalous term 
W"/47T in the central charge will be lost. We will focus only on the logarithmically divergent part of 
the soliton mass renormalization. 
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It can be noticed, that the integral of each term in the last sum in the right-hand 
side of the latter equation over the whole box vanishes identically, since the modes 
in both sectors are normalized to one. However the relevant integral over the fiducial 
interval, discussed in the Sec. [3.9| is not zero, due to different behavior of the modes 
at the edges of the bounding boxQ. In order to calculate the integral over the latter 
interval we again use the quantum-mechanical supersymmetry to express the modes Xn 
in terms of Xm 



X 2 



1/4 \ 2 



UJ1 



Xn (dg-F) Xn 



Xn (p z - F') X 



+ xl 
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XnPP ] y 



X 11 



(6.7) 



Then Eq. ( |6.6p can be rewritten as 



(7~0kink — (^)vac — 9 Z 



W--F 
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(6.8) 



The soliton mass is given by the integral over the fiducial interval, i.e. over z from 



—L/2 + A to +L/2 — A where A is chosen in such a way that L ^> A ^> m 



sec 



the discussion in Sec. |3.9| ). Then it is not difficult to see that the right-hand side in 
Eq.( |6.8| ) vanishes identically within the fiducial interval, and the soliton mass is given 
by the following formula 



M 



W--F' V 



Xn 
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2uj r , 
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2U r , 



(6.9) 



-L/2+A 



This result is perfectly analogous to what we had in the SSG model, except that the 
expression for the "central charge" changes, the one-loop correction to the soliton mass 
is given by 

AM = A(WV> , (6.10) 

with the coefficient of (x 2 ) twice larger than it would be in a calculation of the correction 
to (W) alone. (Note though that the classical part of M is determined by AW.) As a 
result, 

AM 2 ' 



M 



mv 



m 

— In 

4:71 



(6.11) 



This properly reproduces the part of the one-loop correction to the kink mass log- 
arithmic in M r in the non-supersymmetric sine-Gordon model. The finite part is lost 
in the above reasoning because of our formal manipulations with divergent sums. This 
can be readily fixed by considering explicitly an ultraviolet regularization. The finite 



6 A similar behavior can be seen from Fig. |T|: the integrals of the two curves over the whole box are 
equal, while the integrals over the fiducial interval arc clearly different. 
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contribution to the kink mass in fact is well known (see e.g. in |[29|| ) and is equal to 
— m/(27r). We do not present details here, since the goal of this section is to illustrate 
the implementation of the supersymmetric (in the sense of Witten's quantum mechan- 
ics) boundary conditions and the separation of the kink energy from the energy of the 
boundaries. 



7 Anomaly and soliton mass in the SSG model at 
two loops 



We have shown in Sec. [T6] that the anomaly relation ( |3.17| ) is not corrected beyond one 



loop. In the framework of the SSG model we now present two arguments confirming 
this statement. 

The first is based on the relation for the superpotential in the SSG model 

(W"(0)> = ~ <W(0)> , (7.1) 
v 

which holds at all orders of perturbation theory. Using this relation in the expression 
for the anomaly, one finds, assuming that the one-loop expression for the anomaly is 
exact, that at 

2 1 

4:71 

the central charge and also vanish and, thus, the theory becomes conformally invari- 
ant, which is the known behavior of the SSG model at this value of the coupling (see 



e.g. H) 



The second confirmation comes from an explicit two-loop calculation of the kink 



mass, presented below, and its comparison with the known result [32] found under the 
assumption of factorizability of the 5* matrix in this theory. The calculation of the 
kink mass in the SSG model at the multi-loop level was most recently considered in 



Ref. [jHf , where an actual two-loop calculation was also done in detail. We believe that 
our approach, based on the anomalous relation (|3.17|) , provides us with a much simpler 
way of calculating the soliton mass, at least at the two-loop level. The calculation is 
presented in this section. 

According to Eq. (|3.17| ) the mass of the soliton in the SSG model is given by 



M = 2 ( m) + = 2 (l - ^ (W(0)>„ /2 , (7.2) 



where the subscript irv/2 denotes that the quantum average of the superpotential is 
taken over the "flat" vacuum at ((ft) = ttv/2. The relation (|7.1| ) is also taken into account 
in Eq. (|7.2j). Although the calculation of the quantum average of the superpotential 
function in two (or higher) loops is quite straightforward in any given regularization 
scheme in terms of the regulator parameter M r and the bare mass m, we prefer not 
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to present it explicitly, since this part is completely absorbed in the renormalization of 
the particles' mass. Technically the reasoning is most straightforward for the on-shell 
fermion mass, which at the two loop level can be written as 

m F = - (W"(0)) W2 + 5 2 m = M + 8 2 m , (7.3) 

where the equations (|7.2|) and ( |7. 1| ) are used in the last equality, and 82m is the non- 
tadpole contribution to the fermion self-energy given by a single graph of Fig. [| The 



Figure 2: Non-tadpole two loop graph in the fermion self-energy. The heavy line denotes 
the fermion, and the thin lines denote bosons. 



tadpole part of the mass, given by (W" (</>)) has no dispersion, while the finite graph of 
Fig. |2| introduces dispersion of the self-energy in fermion momentum p. For the on-shell 
mass of the fermion, we calculate this graph at p 1 = m and find: 

* 2 ' d 2 k d 2 q 
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(7.4) 

Thus, up to the two-loop terms, inclusively, (i.e. of the relative order v 4 ) the mass 
of the kink is related to the renormalized on-shell mass m of the particles in the SSG 
model as 



M two _i oop = 2 m [v 
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48^2 



(7.5) 



This relation is in perfect agreement with the expansion, up to the appropriate order 
(7 s ), of the known exact relation 
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found [32] within the assumption of a factorized S matrix in the SSG model. 



8 Extra superfields and extended supersymmetry 

In this section we will consider the M = 2 model in two dimensions, obtained by dimen- 



sional reduction of the four-dimensional Wess-Zumino model [191, and its deformations 
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which break M = 2 down to M = 1. A general classification of TV = 2 theories in two 



dimensions is given in [33 



In four dimensions the Wess-Zumino model realizes minimal, M = 1, supersymme- 
try with four fermionic generators. In two dimensions it means that an extra pair of 
fermionic generators is added to superalgebra ( |3.7| ). Let us briefly recall this construc- 
tion putting emphasis on those elements that are specific for M = 2. 

The Wess-Zumino model in four dimensions describes (self-)interactions of one com- 
plex scalar field <£>(x) and one complex Weyl spinor ip a (x) , a = 1, 2 . Both the boson 
and fermion degrees of freedom are united in one chiral superfield []. 

$(x L , 9) = <p{x L ) + V26 a Mx L ) + 2 f(x L ) , (8.1) 

where 9 is a two-component complex Grassmann spinor. The chiral coordinate xl is 
defined as 

so that the supertransformation takes the form 

a -> a + e a , (x L ) a p -> (^l) q/ 3 - 4z , (8.3) 

where a , d = 1, 2 are the Lorentz indices. The Lagrangian of the Wess-Zumino model 
in terms of superfields is 

Avz = \j dA0 ^ + {lj d 2 6W wz (®) + H.c.} , (8.4) 



where Wwz is a superpotential. We choose it to be of the form 

"m 2 _ 2A 



W W z($) = 



$ i:$ 3 

4A 3 



5.5) 



where m and A are real parameters. 

In components the Lagrangian takes the form 

Avz = Ptp+dtf + r*o aa r + f + f + [fw^M - ^w^Mr^a + h. c .} , (8.6) 

where the auxiliary field / can be eliminated by virtue of the classical equation of 
motion, /+ = -W^ z . 

Now, to perform the reduction to two dimensions it is sufficient to assume that all 
fields are independent of x and y, they depend on t and z only. In two dimensions the 
Lagrangian can be presented in the following form: 

1 f - dW d 2 W - 1 



7 For our conventions and notation see Appendix A in Ref. |16 
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We introduced real fields ipi, ipi and fi, where % = 1,2. Up to normalization they are 
just real and imaginary parts of the original fields, e.g, ip = + if 2)/ V2- Summation 
over i is implied. The superpotential W = 2 ImWwz is a function of two variables <pi 
and ip 2 , 

W(^,^ 2 ) = — <pi--<pl + \<pi<f% ■ (8.8) 

The Lagrangian can also be presented in terms of M = 1 superfields $j (see Eq. ( |2.3|) 
for the definition), namely 



C = tjd 2 6 l^DaQiD^i + W($i, $2)} . 



(8.9) 



The presence of extended supersymmetry is reflected in the harmonicity of this super- 
potential, 

- - = for AT = 2. (8.10) 

dcpidfi 

We will consider a more general case of nonharmonic W ((fx, ipz), 

3 



W (<^i, ^ 2 ) = -^r fx - - <Pi + q A y?i0 2 + — <^ 2 



(s.ir 



where p and q are dimensionless parameters. For g ^ 1, p ^ the extended TV = 2 
supersymmetry is explicitly broken down to A/" = 1. 

The classical solution for the kink is the same as in the SPM model with one super- 
field considered earlier; the second field ip 2 vanishes, 



Tfl TfLZ 

<pi = <t>o(z) = ^ tanh_ ^r ' ^2 = 

The topological current £ M now takes the form, 

1 d 2 W ' 



47T C^C^ 



(8-12) 



(8.13) 



For the harmonic superpotentials, i.e. M = 2, the anomaly vanishes. 

The quadratic expansion of the energy density Tt around the solution ( |B.12j) looks 
as follows (compare with Eq. ( |3.45p ): 
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where the following notation is used: 



Xi = <Pi 



X2 = f2 
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(8.14) 



(8.15) 
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The part of this expression containing the fields marked by the subscript 1 is exactly 
the same as in Eq. ( |3.45|) ; the fields with the subscript 2 are new. 

What changes do these new fields introduce? First of all, there is a new fermionic 
zero mode in the field £2, 

£zero = exp(-pmz) ^ 
cosh 29 (mz/2) 

At p = 0, q = 1 it has the same functional form as the old fermionic mode in 77 1. This 
is not surprising because M = 2 supersymmetry places the kink in the supermultiplet 
which contains two bosonic and two fermionic solitons (the second bosonic kink corre- 
sponds to both fermionic modes filled in). What is more intriguing is that the same 
degeneracy remains after M = 2 breaking. If p 7^ and/or q ^ I the M = 2 SUSY 
is replaced by M = 1. Still, two fermion zero modes are present. The Jackiw-Rebbi 
index [R1J] provides a general argument proving the existence of the zero mode in £ 2 - 
Thus, we have two degenerate fermionic solitons, and two bosonic. The degeneracy is 
true only for the ground solitonic states. The M = 2 degeneracy of the nonzero modes 
is lifted at p 7^ and/or q 7^ 1. 

For the nonzero modes of X2, £2 and 772 two new differential operators come into 
play, 

- (d z - 2q\(f) -pm)(d z + 2q\(p + pm)v n {z) = Co 2 n v n {z) , 

-(d z + 2q\(fio+pm)(d z - 2q\(j) - pm)v n (z) = uj 2 n v n (z) , (8.17) 



Again, at p — 0, q — 1 they coincide with the operators of Eqs. fl3.46| ) and ( |3.49|) , due 



to M = 2. This degeneracy is lifted in the general case, but the degeneracy of the 



eigenvalues of the operators (8.17) stays intact. Their eigenmodes are related, 



v n (z) = —(d z + 2q\(f> +pm)v n (z) (8.18) 
Expanding £ 2 in v n (z), X2, V2 in v n (z) and quantizing the coefficients in the same way 



as in Sec. 3.7 we observe the explicit fermion-boson cancellation in the expectation 



value over the soliton state, 

'sol\H(x) -C°(x)|sol) = 0. (8.19) 



The general consideration in Sec. |3.7| proves this to all orders. It is important that in 
proving the fermion-boson cancellation no integration by parts is made. Thus, the state- 
ment of cancellation is insensitive to the particular choice of the boundary conditions. 
(As in Sec. |3T7| they must be M = 1 supersymmetric, of course.) 

So far everything looks very similar to the minimal model treated in Sees. |] and 
|3|. The only distinction is the complex structure inherent to the theory with four 
supercharges. This structure is essential, it is just this complexification that leads to 
the nonrenormalization theorem for the superpotential and for the soliton mass. 
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The relation ( |8.19| ) makes finding the one-loop quantum correction a simple job, 
similar to the one performed in Sec. |j. A somewhat novel feature is that the quantum 
fluctuations of \2 at z — > oo and z — > — oo are different, 
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The result for the soliton mass at one loop is 
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The generalization to the case of arbitrary W(y?i,v?2) presents no problem, 

Mln^ + 2 U 
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where 
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fc=l,2 



S7T 



{M = {0 o (^ = -oo)}} , (8.22) 



and {0o(-s)} denotes the kink solution, v ? i(- 2 ) — 0o(-2), i P2( z ) = 0. We see that there is 
no correction in the case of unbroken Af = 2 supersymmetry, J2 ~Wkk — 0, as it should 
be due to nonrenormalization of the superpotential. 
In the case of the soft breaking by the mass term, 



W = W harmonic + - n 02 



(8.23) 



there is no ultraviolet divergence in the correction to the soliton mass. Moreover, the 
anomaly, presented by the term 2 in the square brackets, also goes away. Indeed, these 
contributions are proportional to Wn + W22 = p- and are the same at both boundaries 
z = ±00. 

It means that we can treat the field $2 as the ultraviolet regulator for $1, provided 
the $ 2 mass is much larger than that of $1, i.e. p 3> In this limit we get 
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(8.24) 



^i=</>o(z=+oo) 



In this way we reproduce again the result ( |4.5p for the soliton mass including the 
anomalous part (up to the substitution 2M r — > p). 



9 Discussion 

The residual supersymmetry in the soliton background field entails a mode degeneracy 
0], in much the same way as in the four- dimensional problem of instantons. This 
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degeneracy manifests itself in the complete cancellation of the quantum corrections, 
order by order in the (small) coupling constant. In the soliton problem this cancellation 
takes place in the quantity (sol|if — 2|sol). The soliton mass coincides with the central 
charge to all orders. This property, once established at the classical level, cannot 
disappear at the level of quantum corrections, at least in the weak coupling regime. 

Surprisingly, this question has caused an ongoing debate in the literature. In many 
works, including some recent ones, it is asserted that BPS saturation is lost at the 
one-loop level. Our conclusions are categorically against this point of view. The reason 
explaining the discrepancy between the (independently calculated in the literature) 
soliton mass and central charge is the omission of the anomalous term in Z. We have 
established the existence of the anomaly (|3.16|) by four independent methods. Although 
our explicit calculations were carried out in one loop, we argue that Eq. ( |3.16| ) is valid 
in the operator form to all orders. In both models considered, SPM and SSG, the 
contribution of the anomaly in the central charge is negative, it is finite (i.e. contains 
no ultraviolet logarithms) and makes Z smaller. Without the anomalous contribution, 
the would-be value of Z apparently exceeds M so \ (e.g. |]l4jl), in sharp contradiction with 
the general principles. In fact, Ref. |14| carries a special section intended to reconcile the 
apparent contradiction between the general requirement that M so \ > Z and the concrete 
result of [TJ]] according to which M so \ < Z. Needless to say that the contradiction is 
automatically eliminated once the anomaly is taken into account. 

To the best of our knowledge, the anomaly in the central charge in two-dimensional 
models with minimal supersymmetry was not discussed previously. However, indirectly 
its existence could have been inferred from well-known results. For instance, the in- 
spection of the exact S matrix in the SSG model (e.g. p^j ) shows that the parameter 
v 2 in the superpotential (see Eq. ( |2.9|) ) and an analogous parameter in the on-shell 
scattering amplitudes are not identical, the latter is different from the former by a fi- 
nite renormalization. In our language this finite renormalization can be traced back to 
the anomaly. 

Much of the controversy in the literature is focussed around the boundary conditions. 
It is not uncommon to choose them in a non-supersymmetric manner, e.g. [[^|13],[T4|]. It 
is quite obvious that Eq. (17) in [14] is incompatible with the residual supersymmetry, 
see Eq. Q2.20| ). Then the manifest cancellation of the boson and fermion corrections is 
lost, and one has to calculate separately the boson and fermion correction. We hasten 
to add that such an approach a priori is fully legitimate, there is nothing wrong with it. 



Moreover, the value of the one- loop soliton mass obtained in this way in Refs. [|I4] ,[1"5H 
and some previous works is perfectly correct. The non-SUSY approach to the kink 
mass calculation has its advantages since the calculation is based on the dynamics 
in the bulk, and the problem of untangling the "boundary condition contamination" 
essentially does not arise. On the other hand, it totally misses all advantages provided 
by supersymmetry. 

The approach we follow, which exploits in full the advantages of supersymmetry, is 
very economic. Nothing comes for free, however. In the SUSY-based approach there 
is a subtle point. If one considers just the soliton mass per se, i.e. the total energy 
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of the system in the box, the boundary conditions may have a drastic impact on the 
result. They distort the soliton state under consideration in the domain adjacent to 
the boundary, see Sec. PT0| . For instance, if the boundary conditions on ^ at z = ±L/2 
are periodic, they force (x(±L/2)) to be the same while in the undistorted soliton 
(x{L/2)) and L/2)) have different signs. This distortion happens even with the 
super symmetric boundary conditions. For instance, in |12| the following condition on 
the derivative of x was suggested: (d z + F')x(L/2) = —(d z + F')x(—L/2), along with 
the antiperiodic boundary condition on x, compatible with ( |2.20| ). In this case the 
boundary conditions enforce {^iP)l/2 = (V ; V ; )--L/2, while for the undistorted soliton 
{4>4')l/2 must be equal to — (^x/j)_l/ 2 . 

It is easy to see that the contamination of the soliton near the boundaries adds to its 
mass a quantity of the same order of magnitude as the one-loop correction. Therefore, 
it is not so surprising that a spectrum of predictions for the corrected soliton mass is 
obtained. These predictions refer to the contaminated soliton, rather than the genuine 
ground-state soliton. It is even conceivable to design a state whose mass at one loop 
is exactly equal to the classical soliton mass. To this end one requires x an d ^2 to 
vanish at the boundary. This is a possible boundary condition, consistent with ( j2.20| ). 
This sensitivity to the boundary conditions is due to the fact that the soliton mass is 
expressed, via SUSY, in terms of the central charge, a quantity defined at the boundary. 

It is clear that the problem with the near-boundary domain is totally unphysical. 
The effects due to the finite-size box must be cut off and discarded. The physical 
ground-state soliton is certainly insensitive to the conditions we impose infinitely far 
away from its center. To solve this problem it is advantageous to consider the local 
form of BPS saturation condition, see Eq. ( |3.42j ). Then, by inspecting the (one- loop 
corrected) soliton profile, or the energy distribution, we can easily figure out which 
contribution to the mass comes from the undistorted soliton and which is related to 
distortions near the boundary. The latter must not be included. This approach is 
close to that suggested in |10[ long ago, the work which was in essence ignored in the 
later publications. Conceptually, we do the same; certainly, we do include the anomaly 
omitted in Ref. |10|| . We expand this approach in various directions, e.g. we consider 
the soliton profile in the minimal (one-superfield) model, and discuss a two-superfield 
model presenting a deformation of the extended M = 2 model. 

If one is not interested in the soliton profile, one can calculate the soliton mass per 
se from the central charge, by shifting the calculation of the latter a little bit, from the 
boundary of the box, inside the fiducial domain. The edge of the fiducial domain is still 
(infinitely) far from the soliton center, and, hence, the central charge is determined by 
the flat vacuum calculation. 

In the M = 2 models the classical expression for the soliton mass stays intact - 
there are no quantum corrections. The corresponding nonrenormalization theorem is in 
essence identical to that valid for the four- dimensional domain walls |17, 18| . The Af = 2 
SUSY can be broken softly down to M = 1 by a mass term of one of the superfields. 
The one-loop correction to the central charge then is finite, since the mass term plays 
the role of the ultraviolet regulator. The expression for the soliton mass is in full accord 
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with the general relation ( |3.16| ). 

One last remark concerning the recent publications [14,15]. To an extent, our study 
was inspired by the thorough analysis of Ref. [14], which clearly exhibited the problem: 
the kink mass received a correction while the central charge stayed intact. In our 
view, the calculations presented in [14] are correct. Our explanation of the results 
obtained in [14] is that they refer to the Hamiltonian corrected by the anomaly and 
to the uncorrected central charge. These quantities, however, are not the ones related 
by supersymmetry. The authors of the recent paper [15] claim a double success: the 
correct value of the kink mass as well as the BPS saturation. Their calculation of the 
kink mass is quite instructive. In particular, it indirectly hints at the anomaly in the 
form of the relation k tan(5g — 8p) = Tn between the boson (8b) arid fermion (8p) phase 
shifts. As for the central charge, it was fixed in [15] to be equal to the kink mass 
by supersymmetry, the assertion ascending to Ref. [4]. The topological nature of the 
charge was not considered. In this way the authors avoided addressing the problem of 
the difference between the uncorrected and corrected forms of the central charge. It 
was implied that renormalization fixes the central charge uniquely - this is the point 
we cannot accept. 

The BPS saturation, as we have demonstrated, is a direct consequence of the surviv- 
ing supersymmetry. It cannot be broken as long as this symmetry is maintained. The 
presence of the anomalous term in the central charge fully resolves the long-standing 
riddle. 
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